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13.  ABSTRACT 


The  contrast  and  size  of  the  target  limit  detection  or  lock-on  range.  Models 
for  the  prediction  of  detection  range  are  evaluated  in  this  report.  The  contrast 
available  at  the  eye  or  the  electro-optical  sensor  is  assumed  to  be  the  limiting 
factor.  The  atmosphere  provides  a  transmission  factor  for  the  target-Uackground 
contrast.  A  model  for  the  prediction  of  contrast  transmission  proposed  by  Duntley 
in  1948  is  examined.  Duntley’s  work  provides  an  analytic  solution  based  on  an 
equation  which  relates  visibility  to  air-transmittance  along  an  inclined  path,  and 
a  table  of  sky-ground  ratios.  An  improved  relation  for  visibility  Is  developed  by 
the  author.  More  realistic  sky-ground  ratios  are  obtained  from  calculated  data  and 
from  flight  test~.< 

The  RRA  Monte  Carlo  model  and  the  AWS  model,  which  predict  contrast  transmission, 
are  compared  for  accuracy.  The  AWS  model  is  also  compared  to  recent  flight  data. 

The  AWS  model  is  shown  to  predict  generally  higher  results  than  the  RRA  data,  due  to 
the  approximations  used.  The  AWS  model  does  provide  a  useful,  fast  tool  for  predic¬ 
tion  of  contrast  transmission. 

A  more  recent  concept  developed  by  Duntley,  the  directional  path  reflectance  R*, 
Is  used  to  describe  the  atmospheric  effects.  R*  provides  a  useful  single  parameter 
for  evaluating  a  situation  to  determine  approach  angles  where  the  effect  of  haze  is 
a  minimum.  The  directional  background  reflectance  t>Ro  must  be  used  with  R*  to 
predict  contrast  transmission.  Graphs  Illustrating  the  application  are  presented. 
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Preface 

The  research  which  went  into  this  report  was  done  at  the  request 
of  the  Operations  Evaluation  Group,  Assistant  Chief  of  Staff,  Studies 
and  Analysis,  USAF.  The  detection  of  tactical  military  targets  is 
approached  from  the  aspects  of  contrast  and  contrast  transmission  by 
the  atmosphere.  An  applied  approach  is  taken  in  which  the  factors 
affecting  contrast  transmission  are  pointed  out,  methods  of  solving  the 
problem  are  presented  and  the  most  promising  approach  chosen.  The 
attempt  has  been  made  to  utilize  existing  models  and  concepts  and  apply 
them  to  this  problem.  Experimental  techniques  have  also  been  investi¬ 
gated  which  may  be  helpful  in  correlating  the  results  of  model  predictions 
and  flight  test  data.  It  has  become  obvious  to  the  author  that  the  rela¬ 
tion  of  model  predictions  to  flight  test  data  is  still  very  inexact  due 
to  the  inaccuracy  of  visibility  measurements  and  the  many  approximations 
which  must  be  made  to  reasonably  handle  the  radiative  transfer  problem. 
This  indicates  a  continued  and  increased  need  for  simple,  reliable  measure 
ment  techniques  during  flight  tests  to  describe  the  physical  quantities 
which  affect  the  operational  problem.  Only  through  these  efforts  can 
predictive  models  be  improved. 

I  would  like  -to  thank  many  people  whose  efforts  made  this  research 
and  its  results  possible.  In  particular,  Colonel  Ed  Battle  who  sponsored 
the  research  allowed  me  freedom  to  pursue  different  approaches, which  has 
been  invaluable  in  completing  the  work.  I  thank  also  Mr.  Tom  Furness  and 
Mr.  Lee  Task  of  the  AF  Aero-Medical  Research  Laboratories  for  their  help 
on  very  short  notice  during  the  measurement  phase  of  the  research  program. 
Captain  Jim  Mardis  of  AFIT  and  Mr.  Ken  Arnold  of  ASD  Technical  Photograph! 


ii 


,’l^J‘7E'~®?3^5^Ty^-SV3P»W355fe^i£^5g£*^^®3l^ij 

1;K^v'r-:-^.,  „, ,  _ 

GEP/PH/72-4 


Branch  also  were  very  helpful  during  the  field  measurements.  Mr.  Lenny  . 
Crouch  and  Capt.  Neil  McQuage  of  the  Air  Force  Avionics  Lab  provided  much 
information  about  the  Air  Weather  Service  Haze  Model.  We  also  had  many 
interesting  and  informative  discussions  on  approaches  to  the  problem,  for 
which  I  am  grateful,  I  thank  also  Dr.  Robert  Fenn  of  the  Air  Force  Cam¬ 
bridge  Research  Laboratories  for  the  interesting  discussions,  guidance 
and  valuable  information  he  provided.  Dr.  Donn  Shankland,  my  advisor  at 
AFIT,  has  provided  some  valuable  insights  and  ideas  which  have  enabled 
me  to  clarify  the  problem.  My  wife  Andrea  has  provided  much  assistance 
in  preparing  drafts  of  the  report  and  more  important,  kindness  and  under¬ 
standing  during  the  many  hours  of  work. 


Edward  A.  Duff 
Captain,  USAF 
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Abstract 

The  contrast  and  the  size  of  the  target  limit  detection  or  lock-on 
range.  Models  for  the  prediction  of  detection  ranges  are  evaluated  in 
this  report.  The  contrast  available  at  the  eye  or  the  electro-optical 
sensor  is  assumed  to  be  the  limiting  factor.  The  atmosphere  provides  a 
transmission  factor  for  the  target-background  contrast.  A  model  for  the 
prediction  of  contrast  transmission  proposed  by  Duntley  in  1948  is  examined. 
Duntley's  work  provides  an  analytic  solution  based  on  an  equation  which 
relates  visibility  to  air-transmittance  along  an  inclined  path,  and  a 
table  of  sky-ground  ratios.  An  improved  equation  for  relating  air-  \ 

transmittance  to  visibility  is  developed  by  the  author.  More  realistic 
values  of  sky-ground  ratio  were  obtained  from  calculated  data  and  from 
flight  tests. 

The  RRA  Monte  Carlo  model  and  the  AWS  Haze  Model,  which  solve  the  | 

radiative  transfer  problem  and  predict  contrast  transmission,  are  compared 

! 

for  accuracy.  The  AWS  model  is  also  compared  to  recent  flight  data  which 
predicts  contrast  transmission.  The  AWS  model  is  shown  to  predict  generally 
higher  results  than  the  Monte  Carlo  data,  due  to  the  allowance  for  only  two 
scattering  events  in  the  AWS  approximations.  The  AWS  model  does  provide  j 

a  very  fast  arid  useful  tool  for  the  prediction  of  atmospheric  effects  on 
contrast  transmission. 

A  more  recent  concept  developed  by  Duntley,  the  directional  path 
reflectance  R* ,  is  used  to  describe  the  atmospheric  effects.  R*  is  a 
calculated  parameter  based  on  the  path  radiance,  downward  irradiance, 
and  air-transmittance.  It  may  be  generated  from  flight  measurements  or 
from  model  calculations.  R*  provides  a  useful  single  parameter  for 
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evaluating  an  operational  situation  to  determine  approach  angles  where  ' 
the  effects  of  naze  will  be  least  on  contrast  transmission.  The  direc¬ 
tional  background  reflectance  bR0  is  shown  to  be  an  important  parameter 
in  the  prediction  of  contrast  transmission.  bR0  must  be  used  with  the 
R*  data  to  predict  the  contrast  transmission.  Sources  of  measured  bR0 
data  are  referenced.  Graphs  illustrating  the  use  of  the  R*  and  bR0  for 
prediction  are  presented. 

The  inherent  target  contrast  is  needed  with  the  contrast  transmission' 
factor  to  predict  detection  or  lock-on  range.  Measurements  of  inherent 
target  contrast  by  photometric  and  photographic  means  are  compared.  A 
simple  measurement  technique  for  field  use  is  recommended.  Radiometric 
measurements  are  presented  to  illustrate  some  spectral  effects  which 
should  be  considered  in  comparisons  between  the  eye  and  different  electro- 
optical  sensors. 
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ATMOSPHERIC  CONTRAST  TRANSMISSION: 
APPLICATION  TO  THE  VISUAL  DETECTION  AND 
ELECTRO-OPTICAL  LOCK-ON  PROBLEM 


I.  Introduction 

Background 

All  visual  systems  detect  a  target  by  an  interaction  of  two  phe¬ 
nomena:  the  size  of  the  target  and  the  contrast  between  the  target  and 
the  background.  Both  of  these  phenomena  have  been  studied  in  detail, 
particularly  during  and  after  World  War  II.  The  size  dependence  is  cne 
of  resolution  and  thus  a  function  of  whether  the  eye  sees  the  target  with 
the  fovea,  or  high  resolution  area,  or  with  the  peripheral  areas  with 
lower  resolving  power.  The  angular  subtense  required  in  each  of  the 
areas  of  the  eye  has  been  described  analytically  as  a  result  of  experi¬ 
ment.  The  contrast  threshold  of  the  eye  is  dependent  on  total  illumina¬ 
tion  present,  the  minimum  detectable  contrast  decreasing  as  the  total 
illumination  increares  until  a  minimum  threshold  is  reached.  Both  of  the 
phenomena  are  also  variable  between  individuals. 

Electro-optical  systems  operate  on  similar  principles  and  have  the 
same  dependencies.  The  target  must  be  large  enough  to  resolve  and  it  must 
present  sufficient  contrast  with  the  background  to  allow  lock-on  and  track¬ 
ing.  To  be  specific,  an  observer  in  an  aircraft  detects  a  target  on  a 
television  display.  He  points  the  electro-optical  sensor  to  center  the 
target  in  the  display  and  places  a  tracking  gate  over  the  target.  The 
electro' ''3  of  the  system  are  used  to  measure  the  received  signal  dif¬ 
ference  from  the  "target"  and  "background".  If  the  contrast  (electron¬ 
ically,  the  signal  to  noise  ratio  or  the  difference  in  signals  divided  by 
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the  background  signal,  where  the  target  is  assumed  to  contribute  only  a 


small  portion  of  the  signal)  is  sufficient,  the  electronics  will  lock-on 


to  the  target  and  will  track  as  the  aircraft  moves.  The  target  size  and 


angular  subtense  is  a  geometry  and  system  resolution  problem.  It  is  very 


important  to  the  detection,  or  lock-on  problem,  but  not  the  subject  of 


this  report.  In  this  regard,  we  consider  the  target  to  be  large  enough  to 


resolve  and  we  consider  the  contrast  as  the  limiting  factor. 


The  problem  to  be  solved  is  the  determination  of  the  contrast  available 
at  some  point  in  the  atmosphere  from  readily  available  meteorological  data 


such  as  surface  visibility  and  cloud  ceiling. 


The  contrast  problem  may  be  divided  into  three  parts.  First,  we  must 


consider  the  contrast  which  exists  between  ihe  target  and  the  background. 


in  the  immediate  vicinity  of  the  target,  the  inherent  contrast.  Seconds 


we  must  consider  how  the  contrast  is  affected  by  the  atmosphere  it  passes 


through  in  getting  to  the  sensor,  the  contrast  transmission  or  atmospheric 


transfer  function.  Third,  we  must  consider  the  required  contrast  at  the 


sensor  for  detection,  or  lock-on  or  tracking,  the  apparent  contrast  thresh¬ 


old.  The  inherent  contrast  may  be  calculated  from  the  reflectivities  of 


the  target  and  the  background,  although  these  are  not  constant  with  wave¬ 


length,  illumination  angle,  and  viewing  angle.  The  apparent  contrast 


threshold  can  be  measured  as  a  function  of  illumination  level  and  specified 


for  the  system.  The  effects  of  the  atmosphere  are  the  most  difficult  to 


solve. 


Work  was  done  by  Koschneider  in  1924  to  predict  the  attenuation  of 
corcrast  along  horizontal  paths  of  sight  (Ref  1:33).  Duntley  extended 
this  work  to  downward- and  upward-looking  paths  of  sight  in  1948  (Ref  2:179). 
Under  the  assumptions  which  Duntley  made,  a  closed-form  solution  was  obtained 
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for  the  contrast  transmission  of  the  atmosphere.  System  designers  and  ' 
systems  analysts  have  subsequently  used  this  solution  by  Duntley  to 
predict  the  effects  of  the  atmosphere  on  contrast  transmission.  The  Air 
Stafr,  in  an  attempt  to  predict  the  performance  of  electro-optical  (EO) 
and  visual  systems,  has  used  the  model,  originally  stated  by  Duntley  and 
reiterated  by  Middleton  (Ref  3:64)  and  by  Bailey  and  Mundie  (Ref  4).  In 
attempting  to  correlate  the  results  of  model  predictions  with  lock-on 
ranges  obtained  in  flight  tests  of  one  EO  guided  system  they  found  large 
discrepancie :  between  some  of  the  results.  It  was  unknown  whether  the 
poor  correlation  was  due  to  scarce  data  about  the  test  conditions  or 
whether  the  model  was  invalid.  In  an  attempt  to  answer  these  questions 
it  was  decided  to  look  in  more  detail  at  the  Duntley  model  and  also  the 
more  recent  models  which  relate  the  effects  of  the  atmosphere  on  contrast 
transmission  to  weather  information. 


Purpose  of  the  Report 


The  purpose  of  this  report  is  to  investigate  contrast  and  contrast 
transmission  by  the  atmosphere  and  to  determine  ways  to  handle  the  pre¬ 
diction  of  lock-on  ranges.  To  do  this  it  was  necessary  to  investigate 
the  Duntley  model  to  determine  the  assumptions  on  which  it  is  based  and 
how  to  determine  the  parameters  which  are  required  to  use  the  model.  Also 
one  had  to  investigate  the  other  models  which  address  the  problem  to  see 
how  they  agree  among  themselves,  and  with  flight  data  which  reports  con¬ 
trast  transmission.  We  then  can  recommend  the  best  model  to  use,  how  to 
use  it,  and  the  limitations  under  which  the  model  must  be  used. 


Scope  of  Report 

The  report  compares  several  models  for  the  prediction  of  contrast 
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transmission  by  the  atmosphere.  It  does  not  attempt  to  >e-sclve  the 

h  * 

•*j>  radiative  transfer  problem.,  but  only  examines  work  which  has  been  done, 

and  applies  it  to  the  problem  at  hand.  The  report  presents  some  measured 
data  on  intrinsic  contrast  and  suggests  how  this  may  be  obtained  simply 
in  future  flight  tests.  It  does  not  examine  the  details  of  any  o.  the 
sensor  systems  involved  or  tneir  requirements  for  apparent  contrast  for 
detection,  lock-on  or  tracking.  In  general,  it  limits  itself  to  the 
photopic  {i.e.  visual)  region  of  the  spectrum  since  the  electro-cptical 
system  is  usually  controlled  by  a  man  who  must  see  the  target  first. 

Several  important  consequences  of  the  different  spectral  responses  of  the 
eye  and  TV  system  are  pointed  out. 

Assumptions 

With  any  model,  assumptions  must  be  made.  These  will  be  pointed  out 
in  the  discussion  of  each  of  the  individual  models.  Along  with  these 
assumptions,  the  apparent  contras  threshold  is  assumed  to  be  a  constant 
for  any  system  to  be  considered,  although  none  of  the  data  here  requires 
that  assumption.  Only  in  the  step  where  lock-on  range  'is  finally  speci¬ 
fied  for  a  given  combination  of  target,  illumination,  and  atmospheric 
condition,  must  apparent  contrast  threshold  be  specified.  Also  it  is 
assumed  that  contrast  is  the  limiting  phenomenon  and  not  target  size  or 
system  resolution.  Again  this  does  not  limit  the  validity  of  the  argu¬ 
ments  presented;  one  must  simply  ask,  when  determining  the  range,  where 
adequate  contrast  is  available  if  tne  target  is  large  enough  to  be  resolved. 

Development  of  the  Report 

In  Chapter  II,  the  basic  equations  of  contrast  and  contrast  trans- 

T*  mission  are  presented  and  important  concepts  are  discussed.  Different 
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approaches  to  the  so'ution;are*laid  out  along  with  some  of  their  assump¬ 
tions.  Duntley's  model  is  examined  in  detail.  In  Chapter  I I I t  an  approach 
for  the  description  of  contrast  transmission  is  chosen  and  discussed  in 
more  detail,  and  the  importance  of  the  directional  background  reflectance 
is  discussed.  Chapter  IV  compares  two  model  predictions  of  contrast  trans¬ 
mission,  and  compares  one  to  measured  flight  data.  Chapter  V  presents 


graphs  which  can  be  used  to  predict  contrast  transmission  under  a  variety 
of  cases.  Some  experimental  measurements  of  intrinsic  contrast  are  pre¬ 
sented  in  Chapter  VI  along  with  important  consequences  of  the  spectral 
nature  of  contrast  and  contrast  transmission.  Chapter  VII  presents  the 
Conclusions  and  Recommendations. 
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II.  The  Problem  of  Contrast  Transmission 
and  Approaches  to  its  Solution 


Radiative  Transfer  Equation 


An  understanding  of  the  physics  of  the  situation  is  helpful  in  under¬ 


standing  the  effects  of  the  atmosphere  on  light.  From  the  conservation 


of  particles  in  phase  space,  one  may  write  the  radiative  transfer  equa¬ 
tion  (Ref  6:27).  The  equation  describes  the  gain  and  loss  of  particles 
from  a  beam  in  direction  n  in  a  small  volume  of  phase  space  resulting 


from  scattering,  absorption,  and  reemission. 


The  equation  may  be  written 


+  Q(r,S2,v,  +  ) 


where 


4>Cr,a,+)  is  the  angular  particle  density,  i.e.,  the  number  of 


particles,  per  unit  volume,  per  unit  solid  angle,  at 


some  position  r,  with  direction  6,  speed  v,  at.  time  t. 


Q(r,fi,v,+)  is  the  source  term  describing  the  number  of  particles  per 


unit  volume,  per  unit  solid  angle,  p«=r  unit  time, 


created  in  the  direction  n,  with  speed  v, 


>(r,Q,n') 


is  the  angular  scattering  phase  function  which  describes 


the  angular  distribution  of  scattered  narticles  at 


point  r.  The  function  is  normalized  such  that 


f  p(r,M')  /  =  f  P  (r.n'.Q) 

J*'  477  ~Ja  477 


dn  =  i 
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8(r,t)  is  the  magnitude  of  the  scattering  coefficient  (units, 

^  per  unit  path  length)  at  position  r  and  time  t  for 

particles  with  speed  v, 

<(r,t)  is  the  magnitude  of  the  attenuation  coefficient  (units, 
per  unit  path  length)  at  point  r  and  time  t  for  par¬ 
ticles  with  speed  v. 

The  first  term  on  the  right  of  the  above  equation  then  describes  the 
net  number  of  particles  which  leave  the  volume  at  r,  per  unit  time,  per 
unit  solid  angle  in  direction  ft,  with  speed  v  and  time  t  which  also  entered 

A 

in  the  same  direction  ft. 

The  second  term  represents  the  particles  per  unit  volume,  per  unit 
time,  lost  from  the  beam  in  all  other  directions,  by  scattering  or  absorp¬ 
tion. 

The  third  term  represents  the  number  of  particles  per  unit  volume, 
per  unit  time,  scattered  into  the  beam  from  all  other  directions  ft' at  time 
r  and  position  r  with  speed  v. 

The  fourth  t°rm  represents  the  total  number  of  particles  per  unit 
volume,  per  unit  time  created  In  the  direction  ft,  with  speed  v,  at  time  t 
and  position  r. 

If  each  of  the  particles  is  a  photon,  the  v  =  c,  the  speed  of  light, 
and  if  we  define 

cftuH|>(r,ft,+)  =  N(r,ft,+)  (3) 

as  the  spectral  radiance,  i.e.  the  energy  per  unit  area,  per  unit  time, 

A 

per  unit  solid  angle  per  unit  energy  interval  at  time  t  in  direction  ft, 
at  position  r  then  we  may  rewrite  Eq  (1)  as 

» 

P' 
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fqr  =  -  cfi*VN~cK(r,+)N(r,n,t) 


4* 


cS(r,t) 

4tt 


,ft',a)N(r,a',+)dQ' 


+  Q'  (r,ft,.Tiio,+) 


(4) 


In  atmospheric  problems,  the  time  rate  of  change  of  the  spectral 
radiance  is  very  slow  compared  to  the  diffusion  time  of  photons  through 
the  atmosphere.  Thus 

U 

Assuming  now  that  the  atmosphere  is  uniform  in  the  x,y  plane  and 
that  variations  occur  only  in  the  z  direction,  we  can  restrict  ourselves 
to  the  case  where  the  transfer  of  radiation  is  independent  of  position 
in  the  x,y  plane.  Then 


where  0  is  the  angle  between  the  path  of  sight  and  the  vertical.  Tnen  the 
transfer  equation  becomes 


|~|cos  e|  —  —  K(Z)N<z,n> 


+ 


B(z) 

4n 


/* 

-/.p(Z,n,Q')N(Z,6')dft'  +  Qf (zfS,tta> 

jq’ 


(6) 


The  physical  meaning  of  the  above  assumptions  is  that  we  neglect 
effects  of  scattered  clouds  which  would  make  the  atmosphere  non-uniform  in 
the  x,y  plane.  Also  we  neglect  any  other  haze  or  fog  effects  which  may  be 
isolated  in  the  atmosihere.  We  also  make  a  flat  earth  assumption. 

We  now  note  that  in  the  first  term  on  the  right  hand  side  of  Eq  (6), 
the  spectral  radiance  N(z,H)  is  made  up  of  two  portions.  It  contains  the 
spectral  radiance  wi.ich  left  some  reference  surface  and  has  not  been 
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attenuated  to  the  point  r  in  the  direction  ft,  and  also  the  spectral 
radiance  which  has  been  scattered  or  emitted  into  the  solid  angle  in  the 
direction  6,  and  has  not  yet  been  attenuated.  We  call  the  spectral 
radiance  at  the  reference  surface  N0(0,ft)  the  inherent  spectral  radiance. 

The  other  portion  of  the  spectral  radiance  we  shall  call  the  path  radi¬ 
ance  N*CZ,ft)  from  the  surface  0  to  the  point  r. 

We  now  consider  our  reference  surface  to  be  the  surface  of  the  earth, 
with  a  coordinate  system  centered  at  the  location  of  our  target  of  interest. 
Then  we  will  be  concerned  with  the  spectral  radiance  Nc(Z+,e,4>)  which  leaves 
the  target  and  background  in  a  small  solid  angle  in  the  direction  ft.  We 
will  be  concerned  with  the  amount  which  reaches  the  position  r,Nr(Z,0,<j>) 
in  that  direction. 

Contrast  and  Contrast  Transmission 

Spectral  contrast  as  defined  here  is  the  ratio  of  the  difference  in 
spectral  radiance  between  the  target  and  background  tc  the  spectral  radi¬ 
ance  of  the  background.  In  this  form  it  is  also  called  universal  contrast. 
Other  forms  of  spectral  contrast  are  also  defined:  contrast  modulation, 
which  is  the  ratio  of  the  difference  in  spectral  radiances  of  target  and 
background  to  the  sum  of  the  spectral  radiances;  and  spectral  contrast 
ratio,  the  ratio  of  target  spectral  radiance  to  background  spectral  radi¬ 
ance.  Only  the  universal  contrast  is  dealt  with  here,  although  Appendix  A 
shows  relations  between  the  different  definitions  of  contrast.  Spectral 
contrast  transmission  is  the  ratio  of  the  spectral  contrast  at  any  point 
along  a  path  of  sight,  the  apparent  contrast,  to  the  spectral  contrast  at 
the  source,  the  inherent  contrast. 

Non-spectral  versions  of  the  above  quantities  are  also  defined  in 
which  the  wavelength  dependence  has  been  integrated  out.  The  derivations 
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which  follow  are  rigorously  true  only  for  monochromatic  light.  The  areas 
where  an  assumption  is  maHe  to  apply  the  results  to  regions  of  the  spectrum 
will  be  pointed  out. 

The  inherent  contrast  depends  on  the  source  and  position  of  illumina¬ 
tion  and  viewing  as  does  the  contrast  transmission.  The  contrast  trans¬ 
mission  also  depends  on  the  general  terrain  albedo  over  which  the  target 
is  viewed.  Albedo  is  defined  here  as  the  ratio  of  the  spectral  reflected 
irradiance  to  spectral  incident  irradiance.  In  addition  contrast  trans¬ 
mission  is  a  function  of  the  aerosol  content  of  the  atmosphere.  Any  other 
absorbing  or  scattering  substances  in  the  atmosphere  will  affect  it.  The 
problem  may  be  laid  out  mathematically  as  follows:  At  any  point  along  a 
path  of  sight  the  spectral  radiance  from  a  target  is  made  up  of  the  image 
forming  light  reflected  from  the  target  and  transmitted  to  that  point. 

No (Z+,6,d>;T(r) ,  and  the  path  radiance  N*(Z,e,<j>),  the  spectral  radiance 
scattered  into  the  path  by  aerosol  and  molecular  components,  and  which 
may  be  called  non-image  forming  light.  Thus, 

Nr(Z,8,4>)  =  N0(Z+,e,$)TCr)  +  N*(Z,6,4>)  (7) 

where  Nr(Z,8,<*)  is  the  spectral  radiance  at  altitude  Z,  receiver  zenith 
angle  8,  and  azimuth  from  the  sun  <j>.  N0(Z+,8,4>)  is  the  spectral  radiance 
at  the  source  r  =  0  into  a  small  solid  angle  about  0 ,<f>  and  T(r)  is  the 
spectral  air- transmittance  to  altitude  Z  along  a  slant  path  r  of  receiver 
zenith  angle  0.  The  geometry  is  illustrated  in  Fig.  1. 

The  equation  holds  for  radiation  from  the  target  or  the  background. 
The  quantities  involved  are  spectral  in  nature  so  that  a  dependence  on 
wavelength  ■'s  implied.  Denoting  the  spectral  radiance  from  the  target  by 
+N  and  from  the  background  by  we  can  compute  the  spectral  contrast  at 
■«*  the  sensor  Cr(Z,8,<}>): 
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cr<z,e,4>>  = 


+Nr(z,e,<f>>  -  bNr(z,e,<i>) 
bNr(z,e,<{.) 


The  spectral  contrast  at  the  source  C0  is 


C0(Z+,e,4>)  = 


+N0 (Z+,0,4>)  -  bN0(z+,e,«) 

bNo(Z+,0,^) 


The  spectral  contrast  transmission  tc  is  the  ratio  of  the  apparent  spec¬ 
tral  contrast  to  the  inherent  spectral  contrast.  It  is  given  by 

C  (z,0,<J>)  ,  , 

Tc(z,0,<j>>  =  c^Tz^mT  ^10i 

Using  Eqs  (8}  and  (9)  and  dropping  the  (Z,0,<j>)  for  brevity, 

tNr  -  bNr  bNo 

T<:  bNr  tNo  “  bN° 

Substituting  from  Eq  (7)  for  the  elements  of  the  first  :erm, 

tNoT  -  bNoT  bNo 

Tc  "  bV  +  N*  '  +No  -  bNo 


Dividing 


bNoT 

Tc  .NT  +  N* 
b  o 


then  dividing  by  the  numerator,  we  get 


This  is  the  basic  equation  for  computing  spectral  contrast  transmission 
through  the  atmosphere.  Rewriting  to  show  the  full  dependence. 


S'*?* 
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Tc(z,e,cf>) 


N*(Z,8  ,<t>> 
bNo(Zt,0,<{>)T(Z>e) 


I 

-f 

5 

I 
:l 

§» 

I 

% 
:4 

'•' 
:■? 

Il 


a 


Here  we  notice  that  the  spectral  contrast  transmission  is  a  function 
of  altitude,  receiver  zenith  angle,  and  azimuth  of  the  path  of  sight  with 
respect  to  the  sun,  and  we  see  the  complex  way  in  which  this  comes  about. 
The  transmission  dependence  is  easily  understood:  as  the  range  increases, 
the  transmission  decreases  due  to  scattering  and  absorption.  The  range 
along  a  slant  path  is  described  by 

R  =  Z  |  sec  6 1  (14) 

where  0  is  the  receiver  zenith  angle  (see  Fig.  1).  This  statement  makes 
the  flat  earth  assumption  and  must  be  corrected  for  long  slant  paths 
(Ref  5:35).  The  spectral  attenuation  coefficient  k(Z)  describes  the 
distance  in  which  a  beam  of  unit  irradiance  is  reduced  to  1/e  of  its 
initial  value  by  scattering  or  absorption  or  both.  The  attenuation  coef¬ 
ficient  k(Z)  =  a(Z)  +  B(Z)  where  ct(Z)  is  the  absorption  coefficient  and 
B(Z)  is  the  scattering  coefficient.  For  a  horizontal  path  k(Z)  =  k„  and 
the  spectral  transmission  is  just 

T(R)  =  exp(  -<„R)  (15) 

Since  the  attentuation  coefficient  is  not  constant  with  altitude,  the 
concept  of  optical  depth  was  evolved.  This  is  defined  by 


WZ)  =  fZ*(Z)  dZ 


The  spectral  air-transmittance  for  any  slant  path  through  the  atmosphere 


is  then  given  by 


T(r)  =  T(Z,0)  =  exp(-x  ,|sec  0|) 
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The  spectral  radiance  from  the  background  is  seen  to  be  a  function  • 
of  target  height,  and  the  angle  of  the  path  of  sight  relative  to  the  sun 
or  other  source  of  illumination.  This  target  height  dependence  enters' 
since  the  amount  of  spectral  irradiance  reaching  the  target  from  the  sun 
or  the  sky  is  a  function  of  altitude.  The  0,<j>  dependence  is  due  to  the 
directional  reflectance  properties  of  all  real  materials.  It  is  only  for 
the  ideal,  Lambertian,  or  diffuse  surface,  that  the  spectral  radiance 
would  be  constant  in  all  directions,  independent  of  the  relation  of  the 
source  of  illumination  and  the  viewing  angle. 

The  dependence  on  altitude  of  the  path  radiance  N*  enters  through  the 
dependence  of  the  number  and  type  of  scatterers,  absorbers  dr  emitters  on 
altitude.  No  single  function  can  describe  this  dependence  although  it  is 
usually  assumed  that  there  is  an  exponential  decrease  in  scatterers  with 
altitude  up  to  some  altitude.  The  dependence  of  N*  on  6  and  <j>  is  due  to 
the  angular  distribution  of  scattering  p(Z,Q,ft' )  with  respect  to  the 
source  of  illumination.  Figure  2  shows  a  Rayleigh  scattering  function  and 
a  representative  aerosol  scattering  function  (Ref  6:21).  The  dependence 


Fig.  2.  Representative  Scattering  Patterns  for 
Isotropic,  Rayleigh,  and  Aerosol  Particles 
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on  wavelength  enters  due  to  the  changing,  with  wavelength,  of  the  scat¬ 
tering  magnitude  (and  shape  for  the  aerosol).  For  the  different  molecular 
(Rayleigh)  components  which  cause  scattering,  this  scattered  intensity  can 
be  written  as 


scattered 


X  **  ( 1  +  cos2$) 


where  X  is  the  wavelength  and  3  is  the  scattering  angle  measured  from  the 
axis  pointed  in  the  direction  away  from  the  source  of  illumination.  For 
aerosol  scattering  the  dependence  is  not  available  in  closed  form  without 
assumptions  about  the  size  and  makeup  of  the  particles.  Since  aerosols 
can  be  smoke  particles,  dust,  water  vapor  haze,  ice  or  snow,  the  specifi¬ 
cation  of  size  and  number  density  is  very  difficult.  It  is  the  task  of 
measurement  to  model  these  parameters  accurately  so  that  calculations  can 

be  made  to  determine  the  scattering  phase  functions.  These  dependencies 

y  "  \ 

listed  above  indicate  the  complex  nature  of  the  path  radiance. 

Following  the  development  by  Duntley  (Ref  7:2-5,  Appendix  A),  the 
last  two  terms  in  Eq  (6)  may  be  tailed  the  path  function  N*(Z,8,<j>). 

Then 


3z[sec  ■9[s-*<Z>N(Z,9,<|>)  +  N*(Z,8,4>)  (19) 

where 

N#<Z,e,4>)  =  /  p(Z,M')N(Z,6')  dQ'  +  Q'(Z,ft'/fi0))  (20) 

47T 

Duntley  then  points  out  that  there  exists  for  each  point  in  the  atmos¬ 
phere  a  unique  spectral  radiance,  the  equilibrium  radiance  Nq(Zj8,<j>),  which 
is  transmitted  unchanged  through  the  small  volume  dz  dfi.  From  Eq  (19),  by 
setting  9N/Oz|sec  0|)  =0, 
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N*(z,e,<^) 

N  ( Z , 6 , cf> >  =  -  (21) 

q  k(Z) 

This  shows  that  at  each  point  along  the  path,  the  spatial  rate  of 
change  of  the  spectral  radiance  is  proportional  to  the  difference  between 
the  spectral  radiance  N(Z,9,<j>)  and  the  equilibrium  radiance  Nq(Z,0,<j>). 
That  is, 

IzjsecHf  =  '  "  Nq(Z,e,<}>)]  (22) 

A  dark  object  will  tend  to  increase  in  spectral  radiance  and  a  light 
object  will  tend  to  decrease,  toward  the  equilibrium  radiance. 


Analytic 


We  now  consider  some  of  the  approaches  to  solving  the  problem,  both 
analytic  and  by  measurement.  We  will  also  consider  the  assumptions  made 
in  each. 


Analytic  Approach  I:  Sky-Ground  Ratio  (Duntley  K  Factor).  Duntley 
originated  the  following  approach  in  1948  (Ref  2:179).  It  has  been  reiter¬ 
ated  and  used,  though  often  incorrectly,  many  times  since  its  origination 
(Ref  3:64  and  Ref  4),  Duntley  develops  the  problem  as  follows  (also  see 
Middleton,  Ref  3:64): 

The  geometry  is  shown  in  Fig.  3a.  In  this  derivation  all  quantities 
are  integrated  over  wavelength  over  the  photcpic  or  visual  region.  To 
avoid  confusion  we  use  the  symbol  B(Z,0,<f>)  to  represent  the  brightness, 
the  photometric  equivalent  of  radiance.  An  equation  for  the  change  in 
radiance  through  the  lamina  dr  is 

=  -  [a(r)  +  3(r)]B(r)  +  B*(r) 


'i  K 


(23) 


If  B*  is  everywhere  proportional  to  8,  the  scattering  coefficient,  then 
Eqs  (23)  and  (24)  hold. 

With  the  above  assumptions  one  may  proceed  to  assume  that 

a(r)  =  a(r)  +  8(r)  =  a0f(r)  (25) 

and 

B*(r)  =  B*(0)f(r)  (26) 

where  a0  is  the  attenuation  coefficient  at  the  surface.  This  assumption, 
as  Middleton  points  out,  restricts  one  to  the  case  whore  absorption  is 
negligible  compared  to  scattering.  This  is  probably  a  good  assumption  if 
large  amounts  of  pollutants  are  not  present  and  also  at  relatively  high 
sun  angles.  We  then  may  write  Eq  (23)  as 

dB(r)  =  [  - a0B( r)  +  B*(0)]f (r)dr  (27) 

Integrating  between  B0  and  BR  and  0  and  R, 
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rBR  _ dR(r) 

*'Bj(  0oB(r)  -  B*(0) 

where  R  is  the  length  of  the  slant  path,  Bq  is  the  radiance  in  the  direc 
tion  of  interest  at  the  source  and  3R  is  the  radiance  at  the  receiver. 
Solving,  we  have 


f 


f (r)dr 


(28) 


I n[a0BR  -  B*(0)]  -  ln[o0B0  -  B*(0)]  =  -a0R 


where 


-  fR 

R  =  ]  f Cr)dr 


which  Duntley  calls  the  "optical  slant  range."  Then 


a0BR  -  B*(0)  =  [a0B0  -  B*(0)]e“a°R  . 

„  B* ( 0 )  aoBe  ~  B*{0)  -aoR 

br  -  “  +  - ^ - 6 

Bp  »  (1  -  e-a°R)  +  B  e~a°R 

R  0O  0 

This  is  seen  to  be  in  the  same  form  as  Eq  (7)  with 


(29) 


(30) 

(31) 

(32) 


B*(0)  ,,  -OoRv 
— - —  (1  -  e  ) 

°o 

equal  to  the  path  radiance  and  e-a°R  equal  to  the  transmittance.  Pro¬ 
ceeding  with  the  same  type  of  development  for  contrast  transmission  we 
obtain 


1  +  jp®  (eao^~T) 
aoDBo 


(33) 


Next  Duntley  goes  back  to  his  differential  Eq  (23)  as  we  did  in  Eq  ('“!)  and 
sets  it  equal  to  zero  to  get  the  equilibrium  radiance: 


‘q 


(r) 


B»(r) 

0(1) 


(34) 
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Hence  if  the  radiance  is  just  equal  to  BfeC«-)/o(r)  it  will  not  be  altered 
with  range.  This  is  the  concept  m?  "optical  equilibrium"  (Ref  2:179). 

By  the  assumptions  made  above  in  Eqs  (25)  a^d  (26), 


B*(r) 


a(r) 


B»(0)  _  R 
t'(0)  ~  ^horizon 


(35) 


since  along  any  horizontal  path  of  sight  B^(r)  is  the  radiance  of  the  hori¬ 


zon  sxy  in  a  particular  direction  provided  chat  the  sky  is  clear. 

What  is  the  direction  along  the  horizon?  Consider  the  downward 
looking  path  cf  sight  as  shown  in  Fig.  1.  The  sun  is  being  scattered  into 
the  smal :  solid  angle  dft  in  the  direction  fi.  From  Fig.  2,  the  angle  is 
described  as  3.  One  must  look  along  the  horizon  in  the  direction  such 
that  the  scattering  angle  is  the  same.  The  sun  must  make  the  same  angle 
with  the  horizon  path  as  it  does  with  the  inclined  path.  The  reason  for 
requiring  clear  skies  is  now  additionally  clear.  This  development  will 
be  valid  only  if  single  scattering  of  the  sun  is  the  prime  cause  of  the 
path  radiance  Bg|^y.  The.i  Eq  (33)  reduces  to 

1 


1  +  (eo°R  .  1) 


(36) 


bBo 


B(<p.  .  =  6)/kB0  is  the  "sky-ground  ratio."  B,  .  must  be  measured 

in  a  particular  direction  which  may  be  specified  as  follows:  If  the 
spherical  coordinate  system  shown  in  Fig.  1  is  set  up,  then  the  cosine  of 
the  scattering  angle  for  a  path  of  sight  in  direction  9,<j>  is  given  by 


(37) 


cos  3  =  cos  0O  jcos  0|  +  sin  0O  sin  0  cos  <j> 

(see  Appendix  B  for  derivation),  where 

0O  is  the  zenith  angle  of  the  source, 

0  is  the  receiver  zenith  angle, 

4>  is  the  azimuth  from  the  sun  <4>  =  0  when  looking  toward  the  sun). 
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and  in  the  direction  along  the  horizon  for  equal  scattering  4>hor|zon  = 


6  =  90°;  thus 


CO?  <{> 


cos  8 


horizon  sin  0, 


(38) 


where  8  is  0°  when  looking  in  the  same  direction  as  the  ongoing  sun  rays. 
Clearly  then,  directions  may  exist  such  that  cos  8  >  sin  90  in  which  case 
cos  4>,  .  >1  and  4>,  .  is  imaginary.  Then  B(<*>.  .  =  8)  cannot 
be  measured.  Figure  3b  shows  a  case  where  the  directions  exist. 

If  the  direction  exists,  and  we  make  all  of  the  above  assumptions, 


then  B(<j> 


hori  zon 


=  8)  can  be  measured  along  with  bB0 (e,4>)  and  the  sky 


ground  ratio  determined.  This  B(4,horjzon  =  8)  is  then  the  non-spectra! 
counterpart  of  the  previously  mentioned  Nq  of  Eq  (21),  the  equilibrium 
radiance  for  a  path  of  sight. 

Duntley  gives  a  table  of  possible  values  of  sky-ground  ratio  under 
some  sky  and  background  conditions.  The  variation  of  sky-ground  ratio 
with  sun  angle  and  path  of  sight  is  not  shown.  As  a  matter  of  fact, 
calling  the  sky-ground  ratio  K,  as  do  Bailey  and  Maudie  (Ref  4)  and  others, 
is  highly  misleading  since  it  implies  that  sky-ground  ratio  is  invariant 
with  source  and  receiver  position.  Actually,  sky-ground  ratio  is  a  func¬ 
tion  of  6  and  (j>,  the  path  of  sight,  80,  the  sun  zenith;  and  the  inherent 
background  reflectance  as  Duntley  points  out  (Ref  2:179). 

This  still  leaves  the  determination  of  R.  Duntley  does  this  by 
assuming  an  “optical  standard  atmosphere"  in  which  the  scattering  and 
absorbing  particle  density  decreases  exponentially  with  altitude.  The 
form  which  Duntley  uses  for  particle  density  is  (Ref  2:187): 

(39) 


N/N0  =  exp  (-  Z/2 1,700) 


where  Z  is  the  height  above  sea  level  in  feet.  Middleton  uses  (Ref  3:74): 
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N/N0  =  exp(-7./30,000)  (40) 

as  do  Bailey  and  Mundie  (Ref  4).  See  Fig.  1  for  geometry.  Since  Z  =  R|cos  0| 
and  assuming  that  the  f(r)  in  Eq  (28)  is  N/N0,  then 

R  =  71,700  Isec  e|^e“RlJC0S  °l/21  ^700  _  e-R2jcos  6j/21 ,700^  (41) 

and  for  the  observer  looking  at  a  target  at  sea  level,  the  expression 


reduces  to 


R  =  21 


,700  [sec  6|^1  - 


-R  |sec  61/21,700^ 


Physically  R  represents  the  slant  path  length  which  has  the  same  number  of 
particles  as  a  path  along  the  surface  of  length  R. 

The  relation  is  then  made  by  Duntley,  using  the  earlier  work  of 
Koschneider  (Ref  1),  between  the  extinction  coefficient  and  the  "meteoro¬ 
logical  range."  If  the  "meteorological  range",  V,  is  defined  as  the  range 
at  which  a  large  black  object  is  just  visible  against  the  horizon  sky,  and 
the  threshold  of  apparent  contrast  for  the  eye  is  accepted  as  .02,  then 
a0  -  3.912/V.  sequent  use  is  made  of  the  relation  a0  =  3/v  in  wh.c'r.  the 
eye  is  assumed  to  have  a  .Qc  threshold  of  contrast.  Then  the  contrast 
transmission  can  be  expressed  in  the  analytic  form 


Tc  " - i^Tf — ~ 

C  1  +  K(e  oK  -  1) 

where  R  is  given  in  Eq  (42)  above,  and 

B(<f>,  .  =  B) 

k  _  horizon 

bB0«,*) 

Bailey  and  Mundie  (Ref  4)  give  a  nice  treatment  showing  some  of  the  effects 
of  different  sky-ground  ratios,  and  meteorological  range  on  contrast  trans¬ 
mission.  The  dependence  of  contrast  transmission  on  air-transmission 
T  =  e_v,C!R  for  "arious  values  of  sky-grcund  ratio  can  be  seen  from  Fig.  4. 
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Note  the  very  strong  dependence  on  sky-ground  ratio  for  a  constant  value' 
of  air- transmittance.  One  key  problem  then  is  to  chooso  the  sky-ground 
ratio,  Duntley  gives  the  data  shown  as  Table  I  (Ref  2:184).  Middleton 
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TABLE  I 

Sky-Ground  Ratios  Suggested  by  Duntley  in  1948 


Sky  Condition 

Ground  Condition 

Sky-Ground  Ratio 

Overcast 

Fresh  snow 

1 

Overcast 

Desert 

7 

Overcast 

Forest 

25 

Clear 

Fresh  snow 

.2 

Clear 

Desert 

1.4 

Clear 

Forest 

5 

observes  that  from  Duntley 's  table  it  appears  that  the  sky-ground  ratio  is 
approximately  the  reciprocal  of  the  surface  albedo  under  an  overcast  sky 
and  1/5  of  the  reciprocal  of  the  albedo  under  clear  skies.  Duntley  gives 
no  justification  for  the  numbers  he  presents.  The  use  of  this  table  and 
Middleton's  observations  on  it  are  a  major  cause  of  error  in  predictions 
using  this  theory.  i 

Calculations  have  been  made  from  the  RRA  Monte  Carlo  data  to  yield  \ 

contrast  transmission  (Ref  8).  By  using  the  computed  contrast  trans-  jj 

mission  and  known  air-transmittance  for  that  model  atmosphere,  the  sky- 

ground  ratio  may  be  computed  from  Eq  (43)  for  different  sun  positions,  '  j 

|  i 

receiver  positions  and  surface  albedos.  The  results  are  presented  graphi-  j  j 

cally  in  Figs.  5  and  6.  In  Fig.  5  the  computed  average  value  over  altitude  j| 

of  sky-ground  ratio  is  plotted  as  a  function  of  albedo  for  three  receiver  |J 
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zenith  angles  and  two  visibilities.  An  isotropic  source  at  the  top  of  • 
the  atmosphere  is  assumed.  The  dependence  which  Duntley  shews  with 
albedo  is  clear  but  the  computed  valuesfor  the  low  albedo  case  of  the 
forest  (albedo  =  .05)  or  intermediate  albedo  of  sand  are  considerably 
lower  than  Duntley  predicts.  Note  also  the  variation  of  sky-ground  ratio 
with  the  amount  of  haze  present.  As  visibility  decreases, the  value  of 
sky-ground  ratio  increases  for  a  constant  albedo.  For  a  plane  parallel 
source  (i.e.  the  sun)  there  is  also  a  marked  variation  with  solar  zenith 
angle  as  shown  in  Fig.  6.  Here  the  computed  average  (over  altitude  and 
azimuth)  sky-ground  ratio  is  plotted  against  solar  zenith  angle  for  three 
receiver  zenith  angles.  Visibilities  of  3  and  25  KM  are  used  and  albedos 
of  .1  and  .9  are  shown.  Note  for  the  low  albedo  the  more  pronounced  effect 
of  solar  zenith.  As  the  sun  approaches  the  horizon,  the  sky-ground  ratio 
increases  sharply.  A  similar  effect  as  in  the  isotropic  case  of  increased 
sky-ground  ratio  with  decreasing  visibility  is  also  noted.  In  general 
the  sky-ground  ratios  are  less  spread  than  Duntley  would  indicate  in 
Table  I.  The  high  albedos  produce  slightly  larger  sky-ground  ratios  than 
predicted  and  the  lower  albedos  produce  lower  sky-ground  ratios.  The 
variation  of  sky-ground  ratio  with  respect  to  azimuth  from  the  sun  can 
also  be  shown  but  the  Monte  Carlo  data  does  not  lend  itself  to  this  type 
of  presentation,  as  will  be  discussed  later. 

Several  other  factors  should  be  pointed  out  about  sky-ground  ratio. 

In  general,  the  sky-ground  ratio  is  not  a  constant  with  altitude,  although 
it  may  be  over  a  limited  altitude  regime.  This  can  be  seen  from  the  RRA 
data,  but  the  altitudes  presented  in  the  original  report  are  .5,  2,  4,  6, 

10  and  50  KM.  These  are  rather  widely  spaced  for  use  at  low  dive  angles 
and  low  altitudes.  Values  of  computed  sky-ground  ratio  from  several  flights 
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Effects  of  Albedo  on  the  Sky-Ground  Ratio  for  an  Overcast  Condition  for  Three  Receiver  Zenith 
Angles.  Data  from  RRA  Monte  Carlo  Calculations. 
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—  3  KM  Visibility 

-  25  KM  Visibility 

X  171°  Receiver  Zenith 
+  130°  Receiver  Zenith 

Z  104°  Receiver  Zenith 
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EFFECT  OF  SOLAR  ZENITH  ANGLE 
AND  ALBEDO  ON  SKY-GROUND  RATIO 

Data  from  Monte  Carlo  Calculation, 
All  Data  Averaged  Over  Azimuth 
and  Altitude 


SOLAR  ZENITH  ANGLE  (deg) 


Fig.  6.  Effects  of  Solar  Zenith  Angle  and  Albedo  on  the  Sky-Ground 
Ratio  for  Two  Visibilities.  Three  Receiver  Zenith  Angles 
Are  Shown.  Data  is  from  Monte  Carlo  Calculation. 
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conducted  in  Germany  in  May  and  June  of  1971  are  shown  in  Figs.  7  through 
10.  The  data  for  these  calculations  were  received  from  Dr.  Robert  Fenn 
of  the  Air  Force  Cambridge  Research  Laboratories  (Ref  9).  The  weather 
and  sun  conditions  under  which  the  flights  were  conducted  are  shown  in 
Table  II. 


TABLE  II 

Geophysical  and  Weather  Conditions  for  Haven  View  Flights 


Flight 

Date 

Solar 

Zenith 

(Avg) 

Altitudes 

Covered 

Sky 

Condition 

A1 bedo 

Visibility 

(Reported) 

c-134 

05/25/71 

55° 

.3-2.4  KM 

Broken 

.07 

33-24  KM 

c-137 

05/28/71 

50° 

.3-3  KM 

Broken 

.04 

15-20  KM 

c-138 

05/29/71 

36° 

.3-5.1  KM 

Scattered 

.03 

13-30  KM 

c-139 

06/03/71 

38° 

.3-5.1  KM 

Scattered 

.04 

30-20  KM 

c-142 

06/06/71 

36° 

.3-5.1  KM 

Scattered 

.05 

7-10  KM 

The  method  of  data  collection  is  described  in  Ref  7  and  will  be  dis¬ 
cussed  briefly  later  in  this  chapter  and  in  Chapter  III.  Given  the  data 
necessary  to  compute  contrast  transmission,  and  the  air- transmittance 
from  in-flight  measurements,  the  sky-ground  ratio  may  be  computed.  The 
values  averaged  over  the  altitude  regime  are  displayed  in  Figs.  7  through 
10.  Figures  7  and  8  show  the  values  of  sky-ground  ratio  in  the  plane  of 
the  sun  plotted  against  receiver  zenith  angle.  Flights  with  similar  solar 
zenith  angles  have  been  grouped.  Figures  9  and  10  show  the  plane  perpen¬ 
dicular  to  the  sun.  There  is  some  agreement  in  shape  between  the  grouped 
data,  but  the  varying  cloud  conditions  could  cause  some  of  the  notable 
differences.  Note  however  the  range  of  values  of  the  sky-ground  ratio, 
again  supporting  the  use  of  lower  values  than  shown  in  Table  I.  Note  also 
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Figure  7.  Sky-ground  ratio  versus  receiver  zenith  angle  in  the  so^ar 
taken  during  Project  Haven  View. 
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that  in  most  of  the  cases  the  sky-ground  ratios  tend  to  approach  a  lower 
value  toward  the' horizon.  In  the  limit  where  Dunley's  theory  reduces  to 
Koschneider's,  the  sky-ground  ratio  would  be  1. 

Despite  the  many  limitations  of  the  sky-ground  ratio  and  tne  fact  that 
it  usually  cannot  be  measured,  as  will  be  shown  in  Chapter  VI,  it  ^an  be 
computed  in  the  way  mentioned  above  from  the  more  complete  models.  Tables 
of  values  generated  from  the  Monte  Carlo  data  for  various  solar  zeniths, 
visibilities,  receiver  positions,  and  albedos,  and  wavelengths  are  shown 
in  Appendix  C.  Sky-ground  ratio  values  computed  from  the  RRA  model  are 
presented  in  Table  III  for  the  following  conditions. 

TABLE  III 


Parameters  for  Sky-Ground  Ratio  Computation  from  Monte  Carlo  Data 


Condition  Visibility 

0 

♦ 

Clear  .55  um  3  KM 

.1 ,  .9  0°,30°,75° 

12  values 

22. 5°, 67. 5° 
112. 5°, 157. 5° 

"  25  KM 

"  0°,30°,70° 

II 

II 

Overcast  .55  um  3  KM 
.65  um 

.1,.3,.6,.S  Isotropic 

II 

- 

"  25  KM 

It  It 

II 

Before  accepting  the  Dunt'iey  model,  with  the  new  sky-ground  ratios, 
some  comments  must  be  made  about  relating  the  visibility  to  the  transmittance 
through  "he  "optical  standard  atmosphere."  The  equation  which  Duntley  gives 
(Eq  42)  and  Middleton  and  Bailey  and  Mundie  repeat,  but  change  21,700  to 
30,000  feet,  appears  to  be  highly  pessimistic  in  predicting  transmission  of 
the  atmosphere.  Work  reported  by  Elterman  (Ref  10:7-1)  uses  an  aerosol 
particle  density  function  of  N/N0  =  exp (Z/ 1.2  KM)  up  to  10  KM.  This  shows 


x 


33 


A* 


6EP/PH/72-4 


...M— 


that  the  aerosol  particle  density  drops  much  more  sharply  than  the  values 
predicted  by  Duntley  in  Eq  (42).  Tins  can  have  a  large  effect  on  the  pre¬ 
diction  of  contrast  transmission  at  the  lower  altitudes.  Elterman  gives 
an  excellent  treatment  of  the  problem  of  transmittance  prediction  as  a 
function  of  surface  meteorological  range  in  Ref  11.  It  appears  that 
Duntley  has  lumped  molecular  (Rayleigh)  scattering  with  aerosol  scattering 
in  coming  up  with  his  optical  slant  range.  Middleton  and  Bailey  and  Mundie 
continue  the  error,  although  they  point  out  that  all  atmospheres  may  not 
look  like  the  standard  atmosphere.  Using  the  development  of  Elterman  (Ref  11) 
and  the  concept  of  optical  thickness  mentioned  earlier,  the  spectral  air- 
transmittance  may  be  developed  in  analytic  form  as  follows: . 

The  transmission  over  a  slant  path  from  Eq  (17)  is 

Thi_h^(Ah,X,Vri)  =  exp  |  ~[Tex+(h2,X,Vn)  -  ,X,Vn>]|sec  0'j  (44) 

where  (h,X,Vn)  refers  to  an  altitude,  wavelength,  and  visibility  condition. 
Tgxt  is  the  optical  thickness  of  the  atmosphere  from  the  surface  to  altitude 


h.  But 


f_,(h,X,Vn)  *  x 


+  T 


+  T. 


'exl h'  lex+D  '  lext  '  lext  , 

R  p  abs 


(45) 


where  R,  p,  and  abs  refer  to  extinction  due  to  Rayleigh  scattering,  aerosol 
scattering  and  absorption  respectively.  By  definition 

•h 


t  ,  (h 
ext 


,X,Vn)  =  J"  6<h,x,vn)dh 


(46) 


for  Rayleigh  particles 


BR(h,X,Vn)  =  eR(h,X)  =  oR(X,h)nR(h) 

since  the  molecular  particle  density  does  not  depend  on  visibility.  oR(X,h) 
is  the  scattering  cross  section  for  molecules.  Also  assuming  uniform  mixing, 
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crR(X,h)  =  oR(X,h=0> 


Text  (h»X>Vn}  =0R(X)J  Vh)Gh  (47) 

R 

Vile  may  now  use  the  U.S.  Standard  Atmosphere,  1962  for  an  expression  for 
number  density  of  molecules  in  the  standard  atmosphere:  N/N0=exp(-Z/9.2  KM). 


Tex+  (h,A)  =  HRBR(X,h=0)[1  -  exp(-Z/HR)]  (48) 

R 

where  HR  =  9.2  KM  is  the  scale  height.  Note  that  this  is  not  a  function  of 
visibility.  Elterman  gives  3R( .55  ym)  =  .0116  1/KM.  Curves  from  Ref  12 
show  the  Rayleigh  aerosol  and  ozone  attenuation  coefficient  as  a  function 
of  wavelength.  The  data  is  presented  in  Fig.  11.  For  aerosols,  a  model 
presented  by  Elterman  (Ref  10:7-3)  shows  a  number  density  for  aerosols 

N/N0  =  e~*-/U2  ^  (49) 

up  to  10  KM  altitude.  The  density  above  that  is  not  easily  approximated. 

The  models  presented  by  Elterman  for  low  visual  ranges  in  Ref  11 
modify  this  slightly,  attributing  the  above  to  clear  conditions  and  add 


haze  below  5  KM.  Here 


N/N,  = 


where  Hp  is  a  scale  height  for  the  particular  visibility  V^,  then 

2(h’x'Y  =  Bp(ho,x,Vn)e-h/HP(vn) 


T  ,  (h,X, V  )  =  H  (V  )B  (h  ,X,V  ) ( 1  -  e' 
extp  n  p  n  p  °  n 


-h/Hp(Vn). 
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WAVELENGTH  (urn) 


culated  Atmospheric  Attenuation  Coefficients  for 
izonl.al  Transmission  at  Sea  Level  in  a  Model 
:ar  Standard  Atmosphere.  Absorption  at  the  Longer 
elengths  Due  to  Hz0  and  CO*  is  not  shown. 
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Neglecting  absorption  and  in  the  photopic  region,  we  may  assume  that  the 
air-transmittance  at  X  =  .55  ym  is  representative  of  the  photopic  region. 
Then 

K0(.55  ym)  =  Bp  +  0R  =  a0  - 
Sc 

%xtp(h'X-V  '  Hf(^r  -  “r)('  -  e-h/HP(V''')  (53) 

For  specific  wavelengths  the  curve  in  Fig.  11  may  be  used  to  find  Bp(X) 
and  pR(X)  or  consult  Ref  10.  Elterman  gives  the  scale  heights  for  the 
visibilities  which  he  models  as  follows: 


Vn  (KM)  Hp  (KM) 


2 

.84 

3 

.90 

4 

.95 

5 

.99 

6 

1.03 

8 

1.10 

10 

1.15 

13 

1.23 

It  appe?-s  appropriate  to  use  l.f  KM  as  the  scale  height  for  visibil¬ 
ities  between  13  and  25  KM  and  also  above  that.  Note  that  the  scale  heights 
for  the  lower  visibilities  are  valid  only  up  to  5  KM  altitude  (Ref  11). 

Above  5  KM  the  decrease  is  no  longer  exponential  (see  Ref  13). 

Up  to  5  KM  in  altitude,  the  transmittance  in  the  optical  standard  atmos- 
phe"e  may  be  written  as 


37 


~«v^»  .^waiiwawK-  • 


GEP/PH/72-4 


•■  T,  .  (Ah,X7V  ) 
hrh2  n 


xp  |-[hrSr(X)(i  -  e-h-'H«)  +  Hp(^p  -  eR(X>) 
-  e'h2/HP(Vn^  -  eR(X)HRfl  -  e"hl/HR) 


-  -¥«)('  -  e-h,/Hp(Vn>)]|sec  6|| 


Combining, 


T.  ,  (Ah,X,V  )  =  exp{  -  BpUm 


h1-h2 


exp  |  -  jeR(X)HR^“hl/HR  -  e“h2/HR^  +  Hp(— 
x  ^“hl/HP(Vr1)  -  e"h2/HP(V^)^|sec  6|| 


P  -  H 


where  h2  <  h2. 


Assuming  that  =  0,  then 


V„(4M.y  =  ®»p  |-[eR<X>HR(l  -  e'h/HR)  +  Hp^2!i  -  eRu>) 
x  (,  -  3-h/HP<V"’)]|sec  e| 


This  expression  may  be  compared  with  Eq  (42)  noting  that  the  Rayleigh  optical 


thickness  is  the  same,  but  the  sizable  correction  of  using  a  scale  height  of 


approximately  1.2  in  place  of  6.41  or  9.2  as  used  by  Middleton  to  describe 


aerosol  extinction. 


Comparison  w.  data  presented  by  Elterman  in  Ref  10  for  a  23  KM  atmos¬ 


phere  and  in  Ref  11  for  a  3  KM  atmosphere  shows  good  agreement  between  the 


formula  given  above  and  the  predictions  based  on  the  data  presented  at  incre¬ 


mental  altitudes.  At  or  below  5  KM  and  at  slant  ranges  less  than  20  KM,  the 


percentage  of  error  between  the  two  methods  is  less  than  14%  even  at  the 
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shallow  dive  angle  of  2.8°.  The  agreement  is  generally  much  better  typi¬ 
cally  less  than  5%  for  both  atmospheres.  Thus  using  the  "modified  optical 
standard  atmosphere"  above  one  can  directly  relate  visibility  to  trans¬ 
mittance  and  use  the  Duntley  sky-ground  ratio  concept. 

Analytic  Approach  II:  Directional  Path  Reflectance  R*.  A  very  inter¬ 
esting  and  promising  approach  to  the  problem  has  been  pointed  out  by  Duntley 
(Ref  7:Appendix  D)  and  also  independently  by  others.  Starting  with  the 
same  basic  equation  for  contrast  transmission,  Eq  (13),  Duntley  suggests 
that  one  divide  the  numerator  and  denominator  of  the  second  term  in  the 
denominator  by  H(Z+,d)T(Z,e) ,  where  H(Z+,d)  is  the  spectral  irradiance 

ir 

on  an  upward  facing,  horizontal  plate  at  the  target  altitude,  d  indicating 
downwelling;  and  Tr(Z,e)  is  the  transmittance  over  the  slant  path  of  in¬ 
terest.  The  equation  now  has  the  form 


[TtN*(Z,e,<|>)/HiZ,d)Tr(Z,e)  ' 
1  +  TTbN0(Z,e,<t>)/H(Z,d) 


-1 


(57) 


where  all  terms  are  spectral.  The  denominator  of  the  second  term,  as 
Duntley  points  out  is  the  "inherent  directional  spectral  reflectance"  of 
the  background  (Ref  7:Appendix  D) .  He  designates  this  by  bR0 ( Z^- ,0,<i>,e  f,«j>  r> 
where  e ',<(>'  refer  to  the  zenith  angle  and  outward  normal  to  the  reflecting 
surface.  Note  that  this  value  would  be  the  same  as  the  albedo  for  a  dif¬ 
fusely  reflecting  or  Lambertian  surface.  Duntley  calls  the  numerator  of 
the  second  term  the  "directional  path  reflectance."  Thus  the  equation  now 
reduces  to  the  much  simpler  appearing  expression. 


Tc(z,e,$>  =  {i  +  CR*(z,e  ,4»)/bRo(z+,0,  ♦)]} 


-i 


(58) 


where 
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Thus  if  one  knows  the  directional  reflectivity  of  the  background,  or  a 
range  of  values  for  it,  then  the  contrast  transmission  is  uniquely  speci¬ 
fied  by  the  directional  path  reflectance.  The  dependence  may  be  illustrated 
by  the  graph  in  Fig.  12  adapted  from  Duntley  (Ref  7: Appendix  D).  Here,  con¬ 
trast  transmission  has  been  plotted  as  a  function  of  directional  path 
reflectance  for  several  values  of  background  reflectance.  Note  that  for  a 
given  background  directional  reflectance,  the  contrast  transmission  decreases 
rapidly  with  directional  path  reflectance  although  it  is  log-linear  over  a 
reasonable  range.  But  note  also  the  relation  of  contrast  transmission  to 
background  reflectance  for  a  fixed  directional  path  reflectance. 

If  one  can  determine  the  directional  path  reflectance,  the  problem  of 
contrast  transmission  may  be  more  simply  treated. 

As  an  added  point  here  it  should  be  pointed  out  that  the  directional 
background  reflectance  is  not  in  most  cases  a  constant  as  for  a  diffuse 


surface,  nor  is  it  purely  specular  like  a  mirror.  It  lies  somewhere  in 
between.  Data  will  be  referenced  that  shows  rather  large  variations  in 
directional  reflectance  of  background  or  target  materials,  thus  showing  a 
strong  variation  of  contrast  transmission  unless  the  directional  path 
reflectance  varies  at  the  same  rate.  This  really  asks  whether  the  path 
radiance  N*(Z,8,<f>)  varies  proportionately  to  the  directional  background 
reflectance.  Duntley  makes  the  statement  that  the  path  radiance  depends 
primarily  upon  the  atmosphere  and  the  position  of  the  sun  and  only  sec¬ 
ondarily  on  the  ground  and  thus  the  directional  path  reflectance  is  only 
mildly  affected  by  the  background  (Ref  7:Appendix  D).  This  appears  rea¬ 
sonable,  since  the  solar  radiance  is  always  much  higher  than  reflected 


„ _ PATH  REFLECTANCE  OR  SEE-FACTOR _ 

Fig.  12.  Dependence  of  Contrast  Transmission  Upon  Directional  Path  Reflectance, 
R*, for  Differing  Values  of  Inherent  Directional  Background  Reflectance 
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radiance  from  terrain,  unless  an  extremely  specular  surface  is  encountered. 
Thus  it  seems  reasonable  that  the  directional  path  reflectance  is  con¬ 
trolled  with  regard  to  the  ground's  effect  by  the  overall  scene  albedo  and 
not  the  directional  background  reflectance.  Thus  the  strong  effect  of 
directional  background  reflectance  on  contrast  transmission  can  be  seen 
from  Eq  (58).  This  approach  will  be  discussed  in  detail  in  Chapter  III. 


6  i; 


Numerical  Calculation  Methods 

The  solution  of  Eq  (13),  the  basic  equation  for  contrast  transmission 
through  the  atmosphere, requires  the  knowledge  of  three  basic  parameters: 
the  path  radiance  N*(Z,e,<(>),  the  transmission  T(Z,e),  and  the  radiance  from 
the  background  bN0(Z,e,<J>).  Duntley  introduced  the  idea  of  the  equilibrium 
radiance  and  optical  standard  atmosphere  to  bypass  the  difficulty  of  mea¬ 
suring  these  values  at  each  point  in  space.  When  the  horizon  direction 
did  not  exist  or  the  skies  were  not  clear,  there  was  no  choice  but  to 
measure  or  do  elaborate  mathematical  calculations.  The  purpose  of  this 
section  will  be  to  briefly  describe  some  of  the  methods  of  calculation 
which  are  used  to  calculate  path  radiance,  transmission,  and  radiance  from 
the  background. 

With  the  advent  of  the  digital  computer,  it  has  become  possible  to 
do  the  complex  calculations  indicated  by  Eq  (6)  or  integral  forms  of  that 
equation.  There  are  several  models  which  exist.  They  have  been  designed 
for  differing  purposes  and  contain  differing  amounts  of  additional  capa¬ 
bilities  above  the  calculations  of  path  radiance,  atmospheric  transmission, 
and  radiance  from  the  target  or  background.  Information  and  in  many  cases 
the  programs  themselves  for  several  of  the  models  can  be  obtained  from  the 
developing  company,  the  Air  Force  Avionics  Laboratory,  Wright-Patterson  AFB, 
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Ohio,  or  from  the  Air  Force  Cambridge  Research  Laboratory,  Hanscom  AFB, 
Massachusetts. 

The  names  of  some  of  the  existing  programs  are  listed  below;  two  of 
the  models  and  the  method  in  which  they  were  used  will  be  discussed  in 
brief  detail. 

1.  Photographic  Reconnaissance  Systems  Analysis  (PRESAC):  Data 
Corporation. 

2.  Evaluation  of  Optical  Haze  on  Atmospheric  Contrast:  Vidya  Cor¬ 
poration  (Ref  14). 

3.  Contrast  Transmission  Data  for  Clear  and  Hazy  Model  Atmospheres: 
Radiation  Research  Associates  (Ref  8). 

4.  MARSAM  Model:  Honeywell  Aeronautical  Division. 

5.  Aerial  Photographic  Energy  Model:  Philco  Ford  Corporation  (Ref  15). 

6.  University  of  Michigan  Atmospheric  Model  (Ref  5). 

7.  Penetration  Survivability  Assessment  Model  (PENSAM):  Honeywell 
Aeronautical  Corporation. 

8.  Atmospheric  Transmission  Model:  Air  Weather  Service  (Ref  16). 

Only  the  third  and  eighth  models  listed  will  be  discussed,  but  infor¬ 
mation  on  all  the  above  may  be  found  in  the  references. 

Contrast  Transmission  for  Clear  and  Hazy  Model  Atmospheres:  Radia¬ 
tion  Research  Associates.  This  model  was  used  to  calculate  the  optical 
parameters  required  by  Monte  Carlo  methods.  The  work  was  done  for  the 
Optical  Physics  Laboratory,  Air  Force  Cambridge  Research  Laboratory. 

Calculations  of  contrast  transmission  are  reported  for  two  model 
atmospheres,  one  with  a  surface  meteorological  range  of  25  KM  (clear), 
and  one  with  a  surface  meteorological  range  of  3  KM  (hazy)  (Ref  8)  A 
diffuse  reflecting  surface  is  assumed  with  albedos  of  .1,  .3,  .6,  and  .9. 
Wavelengths  between  .35  ym,  and  .95  ym  are  used.  Altitudes  of  .5  KM,  2  KM, 

4  KM,  6  KM,  10  KM  and  50  KM  are  reported.  Twelve  downlooking  receiver  angles 


i  a  » Sift*  _ 


ar  =  .*  f".  (a 


GEP/PH/72-4 


and  7  solar  altitudes  are  used.  Azimuths  averaged  over  45°  centered  at 
22.5°,  67.5°,  112.5°,  and  157.5°  are  reported. 

In  addition,  a  plane  isotropic  sour  :  -ne  top  of  the  atmosphere 
simulating  an  overcast  is  assumed  and  Vc  ues  of  contrast  transmission  are 
reported  for  the  same  geometrical  parameters  as  listed  above. 

Table  IV  briefly  describes  the  models  used  for  the  atmospheres,  for 
comparison  with  the  AWS  model  discussed  later.  The  data  has  been  taken 
from  Ref  8. 

Atmospheric  Contrast  Transmission  Calculations: 

The  same  formulation  as  given  above  in  Eqs  (8),  (9),  and  (13)  is  used 
to  define  contrast  and  contrast  transmission.  Path  radiance  N*(Z,8,<|>)  and 
background  radiance  bN0(Z,C  4>)T(Z,8),  (called  "ground  reflectance"  in  Ref  8) 
are  calculated  from  the  light  scattered  intensities,  N(0o,X, albedo, Z, 8, <}>), 
where  the  scattered  radiance  is  the  total  amount  of  spectral  radiance 
arriving  at  Z  in  direction  9S4>.  N(80,X, albedo, Z, 6, 4>)  are  calculated  from 
the  Monte  Carlo  programs  described  in  Ref  17.  N*  and  bN0T  are  calculated 
from  N  as  follows: 

The  scattered  spectral  irradiance  Hs  incident  on  the  background  sur¬ 


face  is 


r2iT  rir/2 

Hs!90,X,a!  bedOjZ^.)  =  J  I  N(90,X, albedo, Z+,0,<J>)  cos  0  sin  8  d0d<}>  (6 

The  background  reflected  spectral  radiance  in  the  direction  n  that  reaches 
the  receiver  at  altitude  Z  is  given  by 

bNo'80,X, albedo, 9>T<Z, 6)  =jHsC0o,X,ai bedo,Zf)  +  exp  Texi.(Z=50  KM)  sec  80j| 

x(alb2&)exp[-  xext(Z>  |sec8|] 


^  / 


(61) 


_ _ 

GEP/PH/72-4 

where  Text(Z--50  KM)  is  the  optical  depth  of  the  total  atmosphere,  (Klbedo/ir) 
is  the  background  reflectance  per  unit  solid  angle  and  exp(-Tex+(Z) jsec  6|) 
is  the  transmittance  from  the  surface  to  altitude  Z.  Then 

N*(6’o,X, albedo, Z,8,4>)  =  N(0o,X,elbedo,Z,6,(i>)  -  bN0 C60 ,X,al bedo,6)T(Z; 9)  (62) 


For  a  plane  isotropic  source  emitting  1  photon/m2  at  an  altitude  of  50  KM, 
the  scattered  intensities  and  .  NoT  were  computed  from  the  Monte  Carlo  data 


for  plane  parallel  source  as  follows: 


N | 5Q^X,a I bedo 


1  ft1'1 

,Z,6)  =  j 

*  7an=f 


N(6o,X, albedo, Z,6,<p)  sin  0o  d0o 


1  P/2 

,N0(X, albedo, Z,0)T(Z, 0)1  =T|  KNo(X,albedo,Z,0)T(Z,0>  sin  0o  d0o  (64) 

D  JISO  L  70o=o  D 


N*S0(X, albedo, Z,0)  =  N!so(X,albedo,Z,0)  -  bN0 (X,al bedo,Z,0)T(Z,0) I  (65) 

I  SO 

Contrast  transmission  may  be  calculated  from  Eq  (13).  Tables  of  contrast 
transmission  are  listed  in  Ref  8,  Vols  I,  II,  and  III. 

A  more  recent  model  has  also  been  developed  by  RRA  for  AFCRL  using  a 
spherical  shell  atmosphere  and  a  "backward"  Monte  Carlo  technique.  The 
author*  4"oicate  that  the  improved  model  is  required  for  twilight  scattering 
studi  sun  is  very  low.  The  model  is  described  in  Ref?  18  and  19. 

AW5  ..Tic  Model.  The  data  summarized  here  for  the  Air  Weather 


Service  model  is. taken  from  "Computer  Simulation  of  Optical  Contrast  Reduc¬ 
tion  Caused  by  Atmospheric  Haze"  (Ref  16),  or  from  the  modified  computer 
program  dated  October  1971,  obtained  from  the  Air  Force  Avionics  Laboratory. 
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Assumptions: 

a.  Horizontally  stratified  atmosphere  (25  layers)  and  composed  of 
atmospheric  gases  ar.d  haze  particles 

b.  Number  density  of  molecules  is  given  by  N/N0  =  exp(-Z/8.430  KM) 

c.  Top  of  atmosphere  is  84.3  KM 

d.  Ozone  absorption  occurs  above  the  top  of  the  atmosphere 

e.  Absorption  of  visible  light  by  other  atmospheric  molecules  and 
haze  particles  is  neglected 

By  neglecting  absorption  in  the  atmosphere  the  extinction  coefficient  at 
any  altitude  is 

Bt(Z)  =  $R'Z)  +  3h(Z)  (66) 

where  T,  R,  and  H  refer  to  total,  Rayleigh  and  haze  scattering  respectively, 
and  the  scattering  coefficient  is 


aT(3„\,Z)  =  or(3,X,Z)  +  aH(B,X,Z>  (67) 

where  B  is  the  scattering  angle  as  discussed  previously.  For  Rayleigh 
scattering 

or(M,Z)  =  *j$r-  o'Z/Hr  (68) 

where  f(B)  is  independent  of  wavelength  and  is  proportional  to  (1  +  cos'1  B) 
aR  is  scaled  such  that  in  the  limit  of  a  Rayleigh  atmosphere  the  surface 
visibility  is  3o6  KM.  For  haze 

[6x(z)]kpx(6) 


°K  = 


4tt 


(69) 


where  p.(8)  is  the  scattering  phase  function  such  that 


/ 

’/4tt 


pa(B) 


dft  =  1 
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the  same  as  the  angular  scattering  phase  function  was  discussed  above  in. 
Chapter  II,  in  cqs  (1)  through  (6). 

Four  atmospheres  are  specified  although  they  may  be  changed  as  required 
in  the  program.  The  data  for  these  atmospheres  is  presented  in  Table  V. 

The  equations  used  in  the  development  of  the  AWS  model  are  given  in 
Ref  16.  The  problem  which  is  solved  is  the  transfer  equation  stated  earlier 
in  Eq  (6).  We  will  describe  here  only  the  basis  of  the  solution. 

As  noted  in  the  list  of  assumptions,  we  have  a  one  dimensional  problem 
with  the  assumption  of  horizontally  stratified  atmosphere.  Then  Eq  (6) 
applies.  Also  k(Z)  reduces  to  6(Z)  with  the  assumption  of  no  absorption 
in  the  atmosphere.  We  also  have  boundaries  at  the  surface  of  the  earth 
i.id  at  the  top  of  the  atmosphere  at  which  we  may  apply  boundary  conditions 
to  solve  our  problem.  The  spectral  irradiance  at  the  "top"  of  the  atmos¬ 
phere  is  well  defined  and  may  be  used  as  an  input.  The  absorption  due  co 
ozone  is  also  obtainable,  sc  this  may  be  taken  into  accr _.,t  above  the  atmos¬ 
phere.  The  spectral  irradiance  reaching  the  surface  of  the  earth  is 
composed  oi  two  portions:  the  direct  solar  spectral  irradiance  which  has 
reached  the  earth  without  scattering  and  the  diffuse  spectral  irradiance 
from  the  sky. 

The  diffuse  irradiance  may  be  thought  of  as  coming  from  single  scat¬ 
tering  events  (where  the  solar  spectral  radiance  is  scattered  once  before 
reaching  the  earth),  and  from  multiple  scattering  events  (which  are  made 
up  of  diffuse  radiance  which  has  been  scattered). 

The  same  ideas  apply  for  any  plane  in  the  atmosphere  but  with  spectral 
irradiances  arriving  from  both  above  and  below. 

The  problem  is  then  solved  in  pieces.  The  direct  solar  irradiance 
reaching  the  earth  is  the  solar  irradiance  after  ozone  absorption  multi¬ 
plied  by  the  transmission  from  the  top  of  t-?  atmosphere  to  the  altitude 
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of  Interest  and  the  cos  0A.  A  diffuse  surface  is  assumed  at  the  earth  • 
so  that  the  reflected  radiance  of  times  the  irradiance  reaching  the 
earth  from  direct  and  diffuse  irradiance  is  reflected  into  a  small  solid 
angle.  This  spectral  radiance  times  the  transmittance  from  the  surface 
to  the  receiver  is  called  the  reflectance  No(Z,0,<j>)T(Z,e)  and  reaches 
the  receiver  at  altitude  Z.  Also  reaching  the  receiver  are  singly  and 
multiply  scattered  radiances  from  the  path  from  surface  to  the  receiver, 
N*(Z,8,<t>). 

The  path  radiance  seen  by  the  observer  is  made  up  of  singly  and  mul¬ 
tiply  scattered  radiances.  The  singly  scattered  radiances  may  be  easily 
handled  using  Eq  (6)  since  we  know  the  spectral  irradiance  from  the  direct 
solar  irradiance  at  any  altitude  along  the  path  from  the  surface  to  the 
observer. 

The  singly  scattered  radiance  reaching  the  earth  which  contributes 
to  the  image  forming  light  may  also  be  handled  from  Eq  (6). 

To  determine  the  multiply  scattered  portions  of  the  radiation  reach¬ 
ing  the  surface  of  the  earth  and  reaching  the  observer,  one  uses  Eq  (6) 
for  an  upward  and  downward  stream  of  diffuse  spectral  irradiance.  Thus 
two  differential  equations  may  be  set  up  and  solved  simultaneously  using 
the  boundary  conditions  to  yield  the  diffuse  spectral  irradiances  moving 
upward  and  downward  at  any  altitude.  The  procedure  related  above  for 
determining  the  singly  scattered  components  is  then  repeated  for  the 
multiply  scattered  components,  using  the  diffuse  irradiances. 

The  program  computes  three  basic  parameters:  the  path  radiance 
N*(Z,0,<j>);  the  air- transmittance  T(Z,6)  from  the  surface  of  the  earth  to 
the  observer;  and  the  inherent  spectral  radiance  N0(Z+,e,(b)  from  a  surface 
of  unit  reflectivity  in  the  direction  ;f  the  observer.  Tnese  three  parameters 
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may  be  used  in  Eq  (13)  by  modifying  the  third  to  coincide  with  the  true  ' 
directional  reflectance  of  the  background.  They  may  also  be  used  to  cal¬ 
culate  the  directional  path  reflectance  R*(Z,0,<f>),  since  the  downwelling 
irradiance  is  just  equal  to  the  spectral  radiance  from  a  unity  reflectivity 
surface  times  tt. 

Measurements:  Ground  and  Flight 

Two  methods  exist  to  estimate  by  measurement  the  effects  of  atmos¬ 
pheres  on  contrast.  The  first  and  direct  way  is  to  measure  the  inherent 
contrast  either  spectral  or  non-spectra!  in  a  given  direction  and  then 
the  apparent  contrast  as  a  function  of  range  along  that  path  of  sight. 

The  ratio  of  apparent  to  inherent  contrast  gives  the  atmospheric  contrast 
transmission.  The  second,  an  indirect  way,  is  to  measure  the  parameters 
to  solve  the  derived  Eqs  (13)  or  (58).  Both  method  '  ere  feasible  but 
involve  different  problems.  For  example,  in  the  direct  method,  one  must 
be  able  to  measure  the  apparent  contrast  between  reasonable  distances  and 
also  be  able  to  describe  the  weather  and  measure  the  visibility  conditions 
so  that  these  measurements  may  be  related  to  other  similar  weather  and 
terrain  conditions.  The  prob1  -  with  the  measurement  then  is  to  maintain 
accuracy  in  the  contrast  measurement  over  a  wide  range,  implying  some 
type  of  zoom  optica!  system  to  maintain  the  same  field  of  view.  Assuming 
».hat  this  could  be  done,  an  extremely  large  number  of  measurements  would 
be  required  to  describe  all  the  possible  approach  angles  and  sun  angles 
and  sky  conditions. 

Contrast  data  measured  photographically  at  200  feet  and  3000  feet  is 
presented  in  Chapter  VI  of  this  report.  A  flight  test  program  to  measure 
apparent  contrast  as  a  function  of  range  from  the  target  under  several 
weather  conditions  is  reported  in  Ref  21.  Two  vidicons  we re  used:  the 
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current  Maverick  vidicon  with  peak  response  in  the  visible  at  .55  ym  and 
a  silicon  vidicon  with  peak  response  in  the  near  IR  at  .75  ym.  Data  was 
reduced  from  the  video  recordings  of  the  individual  TV  lines.  Data  is 
presented  as  apparent  contrast  versus  range  and  an  exponential  curve  is 
fitted  to  the  data.  The  inherent  contrast  and  the  extinction  coefficient 
are  the  free  parameters  for  the  fitting  process.  The  report  should  pro¬ 
vide  a  useful  reference  in  future  attempts  to  reduce  contrast  from  the 
video  recordings. 

The  indirect  approach  is  the  one  presently  being  followed  by  the 
Visibility  Laboratory,  Scripps  Institute  of  Oceanography  in  their  work 
for  Air  Force  Cambridge  Research  Laboratories. 

With  the  concept  of  the  path  function  from  Eq  (20)  and  the  path 
radiance  from  Eq  (7),  Duntley  points  out  that  one  might  experimentally 
determine  the  path  radiance  N*(Z,0,<j>)  by  summing  up  the  products  of  the 
path  function  over  an  incremental  path  Ar  along  the  path  and  the  trans¬ 
mittance  of  the  path  from  that  incremental  path  to  the  sensor  (Ref  7: 
Appendix  A). 


Then, 

m 

N*(Z,8,4>)  =  ]rN*(Z.,e,<j>)T_  (Z,8)Ar 
1  =  1  1 


(70) 


where  N*(Z,6,<}>)  is  the  path  function  and  Tr.(Z,8)  is  the  beam  transmittance 
along  path  r;>  the  distance  from  incremental  path  Ar  to  the  sensor  at  Z. 
N*(Z,8,<j>)  is  given  by 


N#<Z,8,<J>)  =  Nq(Z},8,<j>)8(Zj)  (71) 

where  Nq(Zj,8,4>)  is  the  equilibrium  radiance  for  altitude  Z.  and  6(Z} )  is 
the  total  scattering  coefficient. 
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Nq  is  given  by 


N  <Z,e,4>)  =4 ~J  N(Z{,6\<j>')p(Z,fi',K)dfl' 


(72) 


where  N(Z,er  ,4>f )  is  the  apparent  radiance  of  the  sky  or  ground  for  direc¬ 
tion  6',4>p  and  p(Z,ft,Qf )  is  the  angular  scattering  phase  function  for 
direction  fcr  radiation  entering  in  direction  Q,'  at  altitude  Z{. 

All  the  quantities  here  are  spectrally  dependent.  To  determine  the 
path  radiance  then  for  a  path  of  sight,  one  must  know,  either  experimen¬ 
tally  c:  analytically,  the  radiance  of  the  sky  and  ground  in  all  directions 
from  each  incremental  path  along  the  path  of  sight.  One  must  also  know 
or  measure  the  proportional  directional  scattering  coefficient,  and  the 
transmission  for  each  incremental  path  along  the  path  of  sight. 

The  measurement  techniques  used  to  determine  the  parameters  to  solve 
Eqs  (70),  (71)  and  (72)  are  described  in  Ref  7.  The  measurement  aircraft 
flies  straight  and  level  profiles  at  several  altitudes  and  then  makes  ver¬ 
tical  ascents  and  descents.  During  this  time,  radiometers  with  several 
filters  are  scanning  the  sky  and  ground  to  measure  the  apparent  radiance 
N(Z,e',<jr  )  of  Eq  (72). 

Simultaneously,  a  nephenmeter,  mounted  on  the  end  of  the  wing,  is 
measuring  beam  transmittance  over  an  incremental  oath  for  Eq  (70).  If 
feasible  a  ground  station  is  positioned  along  the  track  cf  flight  to  measure 
the  same  parameters  on  the  ground.  The  data  is  then  interpolated  or  extra¬ 
polated  to  produce  data  at  required  intermediate  intervals.  The  scattering 
functions  used  in  the  cases  reported  here  for  Eq  (72)  are  measured  data 
from  Barteneva  (Ref  7:Appendix  E).  The  path  radiance  may  then  be  calcu¬ 
lated  from  Eqs  (70),  (71)  and  (72).  The  spectral  irradiance  at  the  surface 
may  be  calculated  from 
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H(Z, ,d)  ■*  1  N(0,e',<j>’)cos  0'  dftf  (73) 

•  Zir 

‘he  inherent  directional  background  reflectance  may  then  be  determined 


bN0<zt>M>  ,  , 

bW9'*1  ■  H(Zt,d)/„  <74 

Air-transmittance  is  obtained  by  summing  over  the  incremental  paths. 
Then  contrast  transmission  may  be  determined  from  either  £q  (13)  or 


(58). 

We  note  the  significant  advantage  that  the  indirect  method  is  inde¬ 
pendent  of  any  particular  target. 
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Several  approaches  were  discussed  in  the  previous  chapter.  An 
approach  is  now  chosen  for  solving  the  problem.  It  is  clear  that  the 
measurement  approach  is  highly  desirable  but  very  time  and  resource 
consuming  and  application  to  other  conditions  very  difficult.  Thus 
the  model  approach  is  chosen  here.  The  model  chosen  is  that  in  which 
the  directional  path  reflectance  R*  is  the  key  parameter  rather  than  the 
sky-ground  ratio.  The  latter  still  provides  a  useful  analytic  tool  pro¬ 
vided  the  more  realistic  optical  standard  atmosphere  and  sky-ground  ratios 
presented  in  Chapter  II  are  used.  The  advantage  of  the  directional  path 
reflectance  concept  is  that  it  allov.:  decoupling  of  the  atmosphere  from 
the  directional  background  reflectance.  The  sky-ground  ratio  explicitly 
does  not.  The  directional  path  reflectance  R*  provides  a  single  parameter 
to  describe  the  effect  of  the  atmosphere,  sun,  scene  albedo  combination. 

The  sky-ground  ratio  concept  could  be  adjusted  to  consider  the  individual 
parts,  the  equilibrium  radiance  and  the  background  reflected  radiance,  in 
which  the  same  decoupling  occurs,  but  then  the  formulation  is  more  complex 
since  the  transmittance  must  still  be  considered. 

So  with  the  calculation  of  the  directional  path  reflectance  either 
from  airborne  measured  data  as  described  by  Duntley  et  al .  {Ref  7)  or 
from  a  numerical  model  like  the  RRA  Monte  Carlo  or  AWS  models,  one  can 
use  this  single  parameter  to  investigate  the  ease  or  difficulty  of  con¬ 
trast  transmission  for  any  azimuth  or  dive  angle.  One  can  construct 
surfaces  of  constant  directional  path  reflectance  in  space  and  from 
their  shape  determine  the  best  and  worst  directions  for  contrast  transmission 


S5 


'  . ^ 

...  .*<**#  -*m*tr<Kf 


GtP/PH/72-4 


for  that  atmosphere,  solar  zenith,  and  surface  albedo.  It  must  be  stressed 
that  the  contrast  transmission  is  still  not  uniquely  determined  since  ^he 
directional  background  reflectance  must  be  used.  The  decoupling  mentioned 
above  aids  greatly  in  simplifying  the  problem  of  determining  the  atmospheric 
effects.  Some  of  these  effects  can  be  seen  “rom  Figs.  13  through  18.  Figure 
13  plots  the  directional  path  reflectance  versus  azimuth  from  the  sun  for 
a  clear  atmosphere  (23  KM)  but  with  two  solar  zenith  angles.  A  wavelength 
of  .55  ym  and  an  albedo  of  .06  have  been  used.  Note  the  stronger  effect 
of  azimuth  from  the  sun  as  the  solar  zenith  angle  increases.  The  data 
has  been  calculated  from  the  AWS  program.  Note  also  the  increase  in  direc¬ 
tional  path  reflectance  as  the  zenith  angle  of  the  sun  increases,  implying 
better  seeing  conditions  in  general  for  the  higher  sun.  In  addition,  the 
effect  of  azimuth  appears  strong  for  only  a  limited  azimuth  range,  while 
there  is  a  rather  broad  area  in  w^'ch  the  effect  is  not  too  great.  Now 
going  to  Fig.  14,  the  contours  of  constant  directional  path  reflectance 
have  been  plotted  versus  altitude  and  horizontal  range  for  a  sun  angle  of 
42°  a~d  an  albedo  of  .06.  Figure  15  shows  the  same  type  of  plot  except 
the  solar  zenith  has  been  increased  to  78°.  Both  plots  are  presented 
for  an  azimuth  cf  90°  from  the  sun.  Note  in  general  the  increased  slant 
range  as  the  dive  angle  increases.  This  improvement  is  due  to  an  increase 
in  transmission  with  dive  angle  and  a  decrease  in  path  radiance.  Again 
the  data  has  been  generated  from  the  AWS  program.  The  contours  differ  as 
the  azimuth  is  changed  as  noted  from  Fig.  13. 

Importance  of  Directional  Background  Reflectance 

Once  the  directional  path  reflectance  is  known,  one  must  have  the 
angular  dependence  of  the  directional  background  reflectance  to  determine 
the  contrast  transmission.  Many  times,  since  this  is  not  known,  a  diffuse 
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Contours  of  Constant  Directional  Path  Reflectance  R*  for  a  Solar  Zenith  Angle  o 
Azimuth  of  90°.  Albedo  is  .06  and  X  =  .55  pm.  Values  of  Contrast  Transmission 
Assuming  hRo  *  Albedo.  Data  from  AWS  Model. 
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15.  Contours  of  Constant  Directional  Path  Reflectance  R*  for  a  Solar  Zenith  Angle  of  78°  and  an 
Azimuth  of  90°.  Albedo  is  .06  and  X  is  .55  ym.  Values  of  Contrast  Transmission  are  Shown 
Assuming  hR0  =  Albedo.  Data  from  AWS  Model. 


surface  will  be  assumed,  in  which  case  the  directional  background  reflec¬ 
tance  is  just  the  scene  albedo  used  to  calculate  the  directional  path 
reflectance.  In  the  latter  case,  the  three  figures  (13,  14,  and  15)  de¬ 
scribe  exactly  the  shapes  of  the  contours  of  constant  otrast  trans¬ 
mission.  In  that  case  Eq  (58)  or  Fig.  12  may  be  used  to  convert  the 
contours  of  constant  directional  path  reflectance  to  ones  of  constant 
contrast  transmission.  Similarly,  the  'igures  may  be  rescaled  in  terms 
of  contrast  transmission.  This  is  shown  in  Figs.  14  and  15.  But  if  the 
directional  background  reflectance  is  not  a  constant,  then  the  shapes 
of  the  curves  above  are  changed.  From  the  discussion  in  Chapter  II,  we 
saw  the  strong  effect  of  directional  background  reflectance i  Here  we 
reference  data  collected  by  Duntley  and  others  for  AFCRL  during  flight 
tests  which  shows  the  directional  background  reflectance  variation  with 
look  angle  and  azimuth  from  the  sun  for  a  particular  type  of  terrain 
(Ref  22:795).  The  data  was  taken  using  a  photopic  filter  so  that  it  is 
an  integrated  reflectance.  It  must  be  first  noted  that  this  data  is 
based  on  an  average  reflected  energy  from  a  large  area  of  not  necessarily 
uniform  terrain  and  that  the  value  which  should  be  used  in  Eq  (58)  is 
the  directional  reflectance  of  the  background  in  the  immediate  vicinity 
of  the  target.  But  the  variation  shown  should  be  representative  of 
terrain  of  that  type.  The  data  is  presented  in  Figs.  10  and  17.  It  was 
taken  with  solar  zenith  angle  of  42°  and  the  albedo  was  .06.  If  thir 
data  is  then  used  instead  of  assuming  a  constant  value  of  directs 
background  reflectance,  it  can  be  applied  through  Eq  (58)  or  Fig.  12 
to  generate  the  contours  of  constant  contrast  transmission  shown  in  Fig. 
18.  Here  we  assume  that  the  R*  data  is  representative  of  the  photopic 
or  visible  region.  Note  the  pronounced  change  in  shape  between  Fig.  14 
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and  Fig.  18.  Other  data  measured  by  Duntley  indicates  the  same  type  of  ■ 

M.  i* 

strong  deperdence  with  look  angle  and  azimuth  (see  Table  VI).  The  impor- 
lance  of  knowing  tho  directional  background  reflectance  is  clearly  illus¬ 
trated,  particularly  for  low  dive  angles.  Several  sources  of  this  data 
are  available.  The  sources  listed  here  are  primarily  those  done  by  Duntley 
and  others  for  AFCRL.  The  University  of  Michigan,  Target  Signature  Anal¬ 
ysis  Center  has  compiled  a  data  bank  for  NASA  which  contains  directional 
background  reflectance  data  for  many  background  materials.  In  Table  VI 
t'>e  target  or  background  material,  the  sun  angles  and  the  sky  conditions 
are  listed  alGng  with  the  reference  where  the  data  appears. 
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reasonable  agreement  to  the  Monte  Carlo  model,  comparison  will  be  made  to 
flight  data  collected  during  Project  Haven  View. 

The  more  direct  calculational  approaches  like  the  AWS  model  and  the 
University  of  Michigan  model  possess  superiority  over  the  Monte  Carlo 
method  in  making  rapid  calculations  on  the  computer.  A  rough  estimate  has 
indicated  that  the  fWS  program  runs  360  cases  in  the  time  the  Monte  Carlo 
calculations  program  runs  1  case.  Here  a  case  is  defined  as  a  particular 
visibility,  solar  zenith,  albedo,  wavelength,  altitude,  azimuth  and  re¬ 
ceiver  zenith  angle  position.  This  can  be  very  important  if  one  is 
attempting  to  describe  a  complex  situation  and  trying  to  optimize  a  situa¬ 
tion  for  best  seeing  conditions.  One  could  easily  inciuse  the  i~st  program 
as  a  subroutine  in  an  overall  program,  while  the  Monte  Carlo  calculation 
method  would  almost  be  out  of  the  question. 

T  a  question  then  arises  as  to  how  the  results  of  the  faster  AWS  pro¬ 
gram  can  compare  to  the  results  predicted  by  the  more  exact,  and  by  its 
nature  more  time  consuming,  Monte  Carlo  method.  The  RRA  Monte  Carlo  data 
for  the  calculations  are  from  Contrast  Transmission  Data  for  a  Clear  and 
Hazy  Model  Atmosphere  (Ref  8),  already  described.  The  AWS  data  has  been 
run  crom  the  computer  program  as  described  in  "Computer  Simulation  of  Op¬ 
tical  Contrast  Reduction  Caused  by  Atmospheric  Aerosol,"  (Ref  16)  previously 
cited.  Since  the  computer  program  was  not  available  for  the  Monte  Carlo 
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calculations,  and  the  AWS  program  was,  the  AWS  data  was  generated  to  match 
as  closely  as  possible  the  conditions  already  generated  by  the  Monte  Carlo 
method  and  then  the  comparison  was  made.  The  comparison  could  not  be  exact 
without  some  extensive  although  not  difficult  modification?  to  the  AWS  pro¬ 
gram.  These  were  not  done.  Table  VII  indicates  the  value  of  each  parameter 
used  in  the  Monte  Carlo  calculations  and  the  comparable  parameter  used  in 
the  AWS  program. 

ms  noted  from  Table  VII  the  major  differences  occur  in  the  altitudes 
and  azimuths  used  for  the  calculations  although  some  other  differences 
exist  in  the  scattering  phase  functions  as  may  be  detected  from  the  dis¬ 
cussion  in  Chapter  *T.  The  Monte  Carlo  method  because  of  its  statistical 
nature  requires  intervals  over  which  the  photons  are  collected,  while  the 
AWS  method  allows  for  calculations  at  discrete  points.  Also,  the  AWS  model 
requires  inputs  in  whole  KFT  so  the  altitude  match  was  not  exact.  The 
largest  error  in  altitude  match  is  for  the  l«w  altitude  case  in  which  the 
AWS  data  point  is  20%  higher  in  altitude.  Averaging  over  azimuth  can  also 
have  a  strong  effect.  As  can  be  seen  from  Fig.  13,  already  discussed,  the 
directional  path  reflectance  and  thus  the  contrast  transmission  can  vary 
sharply  in  the  region  between  an  azimuth  of  0°  and  about  45°  and  also  near 
180°  depending  on  solar  altitude.  Thu'  araging  over  the  interval  0°-  45° 
and  135°  -  180°  and  calling  these  the  mid-value  of  angle  for  comparison  may 
lead  to  error. 

In  comparing  the  two  mo>is,  the  following  formula  was  used  to  deter¬ 
mine  percentage  difference: 
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TABLE  VII 


Parameters  Used  for  Monte  Carlo  -  AWS  Comparison 


Parameter 

Monte  Carlo 

Air  Weather  Service 

Visibility 

3  KM 

3  KM 

25  KM 

25  KM 

Aerosol  extinction  modified  to  be 

similar  to  Monte  Carlo  from  Ref  8 

for  both  visibilities 

Wavelength 

.55  uni 

.55  ym 

A1 bedo 

.1 

.1 

.9 

.9 

Solar  zenith 

(P 

0° 

30° 

30° 

70°{25  KM  Vis) 

70°(25  KM  Vis) 

75°(  3  KM  Vis) 

75°(  3  KM  Vis) 

Azimuth 

22. 5° (0-45°) 

?3° 

67. 5° (45-90°) 

68° 

112.5°(90-I35°) 

113° 

157.5°(135-180°) 

158° 

A1 ti tude 

.5  KM 

2  KFT( .61  KM) 

2.0  KM 

7  KFT (2.13  KM) 

4.0  KM 

13  KFT(3.96  KM) 

6.0  KM 

20  KFT (6.10  KM) 

iO.O  KM 

33  KFT (10. 05  KM) 

Receiver  angles 

170.533° 

171° 

164.267° 

164° 

157.667° 

158° 

148.217° 

148° 

138.583° 

138° 

130.550° 

130° 

123.667° 

123° 

116.750° 

116° 

110.383° 

110° 

104.483° 

104° 

All  receiver  angles  are 

actually  the  center  of 
a  small  interval  in  the 
Monte  Carlo  calculation; 
see  Ref  8  for  exact  inter¬ 
val 
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computed  sky-ground  ratios.  Thus  a  negative  percentage  implies  that  the 
AWS  model  calculates  too  high,  and  a  positive  percentage  implies  that  the 
AW5  calculates  too  low  (here  we  make  the  assumption  that  the  '  ante  Carlo 
method  is  more  accurate). 

25  KM  Case 

Let  L_  first  compare  the  air-transmittance  calculation.  Here  the 
aerosol  profiles  used  are  very  similar  as  seen  from  the  previous  discussions. 
Both  are  drawn  from  the  same  scarce,  Elterman  (Refs  10  and  13).  Fhe  AWS 
data  was  adjusted  to  produce  a  surface  visibility  of  25  KM  used  on  the  RRA 
model . 

Comparison  of  the  contrast  transmission  data  is  illustrated  in  Figs. 

19  through  28.  Figure  19  shows  the  sun  at  0°  solar  zenith  (directly  overhead). 
Contrast  transmissions  are  piotted  for  both  models  as  a  function  of  altitude 
for  three  receiver  zenith  angles.  The  RRA  data  averaged  over  azimuth  is 
plotted  since  there  is  large  statistical  variation  with  azimuth  which 
should  not  be  present.  The  results  for  the  two  steeper  angles  show  that 
the  AWS  model  is  more  optimistic  than  the  RRA  model  and  the  percentage 
difference  increases  with  altitude.  At  the  shallow  look  angle,  the  agree¬ 
ment  is  good  although  the  AWS  data  tends  to  be  more  optimistic.  For  an 
albedo  of  .9,  all  other  conditions  being  the  same,  the  data  is  plotted  in 
Fig.  20.  For  ail  receiver  angles,  the  AWS  model  is  more  optimistic;  again 
the  difference  increasing  with  altitude.  For  a  relatively  high  sun,  solar 
zenith  angle  of  30°  and  an  albedo  of  .1,  contrast  transmission  for  the 
two  models  is  plotted  for  an  azimuth  of  113°  i.n  Fig.  21.  Here  the  agree¬ 
ment  is  much  better  for  the  two  shallow  look  angles  but  poor  for  the  steep 
angle.  Thi:  large  difference  appears  to  be  due  to  the  statistical  varia¬ 
tion  in  the  Monte  Carlo  method  rather  than  physical  reasons  since  comparison 
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Comparison  o'  Calculation  of  Contrast  Transmission  by  the  AK3 
and  RRA  Monte  Carlo  Methods.  The  Wavelength  is  .55  ym.  Solar 
Zenith  Angle  is  0°, Albedo  is  .1.  The  Surface  Visibility  is  25  KM 
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25  KM  VISIBILITY 
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Comparison  of  Calculation  of  Contrast  Transmission  by  the  AWS 
and  RRA  Monte  Carlo  Methods.  The  Wavelength  is  .55  pm,  Solar 
Zenith  Angle  is  0°, Albedo  is  .9.  The  Surface  Visibility  is  25  KM. 
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for  look  angles  of  164°  and  158°  shows  a  maximum  of  11%  error  under  these 
same  conditions.  The  general  trend  continues  to  show  an  it. crease  in  dif¬ 
ference  with  altitude  with  AWS  predicting  slightly  higher  values.  In 
Fig.  22  the  effect  of  azimuth  is  examined.  Here  for  an  altitude  of  2  KM 
and  for  two  look  angles,  contrast  transmission  is  plotted  against  azimuth 
for  a  solar  zenith  of  30°  and  an  albedo  of  .1.  The  Monte  Carlo  data  is 
presented  in  bar  graph  fashion  due  to  the  averaging  over  a  45°  azimuth  in 
their  presentation.  The  agreement  in  shape  is  reasonable.  The  12%  dif¬ 
ference  for  the  look  angle  of  131°  and  azimuth  of  23°  appears  to  be  accounted 
for  by  the  statistical  nature  of  the  Monte  Carlo  program,  since  this  con¬ 
figuration  is  near  the  minimum  of  the  scattering  lobe  for  primary  scattering, 
and  path  radiance  should  be  a  minimum  assuming  single  scattering.  Figure  23 
presents  the  same  conditions  as  in  Fig.  21  but  for  an  albedo  of  .9.  Here 
the  AWS  model  is  again  optimistic, and  in  comparison  to  the  .1  albedo  the 
difference  between  the  two  models  is  greater  for  the  higher  albedo.  This 
is  probably  due  to  the  increased  upwelling  radiation  making  multiple  scat¬ 
tering  more  likely.  The  Monte  Carlo  method  should  handle  multiple  scattering 
more  accurately  than  the  AWS  model.  The  Monte  Carlo  data  for  the  look  angle 
of  171°  seems  unrealistically  low,  particularly  since  it  is  well  below  the 
131°  data.  Data  at  1G4°  and  158°  show  a  maximum  of  10%  difference  from  the 
AWS  data.  Figure  24  presents  similar  conditions  to  Fig.  22  but  with  an 
albedo  of  .9.  The  AWS  data  is  still  optimistic,  more  so  than  for  the  albedo 
of  .1  and  probably  due  to  increased  multiple  scattering  mentioned  above. 

The  insensitivity  of  the  AWS  duta  to  azimuth  compared  to  the  Monte  Carlo 
data  may  be  also  due  to  the  Monte  Carlo  superiority  in  accounting  for  mul¬ 
tiple  scattering  events.  Figure  25  shows  the  case  of  solar  zenith  angle 
70°  and  al b<~  .1.  The  agreement  for  all  anales  is  reasonably  good  with 
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Fig.  23.  Comparison  of  Calculation  of  Contrast  Transmission  by  the  AWS 
and  RRA  Monte  Carlo  Methods.  The  Wavelength  is  .55  ym.  Solar 
Zenith  Angle  is  30°, Albedo  is  .9.  The  Surface  Visibility  is  25  KM. 
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Fig.  25.  Comparison  of  Calculation  of  Contrast  Transmission  by  the  AWS 
and  RRA  Monte  Carlo  Methods.  The  Wavelength  is  .55  ym.  Solar 
Zenith  Angle  is  70°, Albedo  is  .1.  The  Surface  Visibility  is  25  KM 
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the  AWS  model  again  predicting  more  optimistically  as  the  altitude  increases. 
Figure  26  shows  the  variation  with  azimuth.  The  agreement  here  is  also  quite 
good,  although  the  Monte  Carlo  point  for  look  angle  of  131°  is  probably  too 
high.  Figures  27  and  28  present  the  same  information  for  an  albedo  of  .5. 
Previous  tre  !s  are  repeated  with  AWS  more  optimistic  and  increasing  slightly 
with  altitude.  The  significant  increase  in  percentage  difference  is  not  as 
notable  here  as  it  was  in  the  comparison  of  albedo  for  the  solar  zenith  of 
30°. 

One  can  then  conclude  that  the  predictions  of  contrast  transmission 
by  the  AWS  model  are  comparable  to  those  of  the  more  time  consuming  and 
statistically  variable  Monte  Carlo  method  for  the  clear  visibility  case 
compared  here.  This  includes  angles  down  to  solar  zeniths  of  70°  and  look 
angles  down  to  104°.  The  AWS  model  tends  to  be  more  optimistic  than  the 
Monte  Carlo  method  and  the  percentage  difference  increases  with  altitude. 

An  increase  in  percentage  difference  with  albedo  was  also  noted  for  the 
high  sun  angles,  but  not  for  the  low  angles. 

3  KM  case 

Due  to  the  difference  in  aerosol  profiles  used  by  the  RRA  Monte  Carlo 
method  and  the  AWS  model,  the  AWS  model  aerosol  extinction  profile  was  changed 
so  that  it  reflected  exactly  the  values  used  in  the  RRA  model,  up  to  and  in¬ 
cluding  3  KM  altitude.  As  noted  previously,  both  models  use  the  same  data 
above  that  altitude. 

Similar  sun  geometries  and  albedos  are  compared  for  the  3  KM  case. 

The  results  are  shown  in  Figs.  29  through  36. 

Figures  29  and  30  show  the  case  for  solar  zenith  angle  of  0°  and 
albedos  of  .1  and  .9  respectively  for  three  receiver  look  angles.  Note  the 
good  agreement  for  the  low  albedo  but  the  large  difference  for  the  high 
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Fig.  27.  Comparison  of  Calculation  of  Contrast  Transmission  by  the  AWS 
and  RRA  Monte  Carlo  Methods.  The  Wavelength  is  .55  ym,  Solar 
Zenith  Angle  is  70°, Albedo  is  .9.  The  Surface  Visibility  is  25  KM. 
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albedo  case.  For  the  low  albedo  the  AWS  model  is  more  pessimistic  than, 
the  RRA  Monte  Carlo  data,  while  for  the  high  albedo  the  reverse  is  true. 
Figures  31  and  32  show  the  cases  for  solar  zenith  of  30°  at  an  azimuth  of 
113°,  and  the  low  and  high  albedos.  Again  the  agreement  is  quite  good  for 
the  low  albedo  case  with  the  Monte  Carlo  generally  slightly  higher  although 
the  data  does  cross  for  the  look  angle  of  171°.  For  the  high  albedo  case 
the  AWS  data  predicts  significantly  higher  results  than  the  Monte  Carlo 
data.  The  azimuth  dependence  is  shown  in  Figs.  33  and  34.  Here  the 
data  is  plotted  for  an  altitude  of  2  KM.  The  agreement  in  shape  is  rea¬ 
sonable,  but  note  the  relative  insensitivity  to  azimuth  of  the  AWS  data. 
This  insensitivity  was  also  noted  in  the  25  KM  da.ta.  Figures  35  and  36 
show  the  case  of  low  sun,  a  solar  zenith  angle  of  75°,  at  113°  azimuth, 
for  the  two  albedos.  Note  that  the  same  trend  occurs  as  for  the  higher 
suns  but  now  the  agreement  is  better  for  the  high  albedo  while  for  the 
low  albedo  the  AWS  data  is  significantly  lower  than  the  Monte  Carlo  data. 
Figures  37  and  38  show  the  azimuth  variation  for  this  sun  angle.  Again 
the  AWS  data  is  quite  insensitive  to  azimuth  change.  This  difference  seems 
to  be  significant  since  in  comparison  to  the  25  KM  da;a  in  Figs  26  and  28 
somewhat  more  change  might  be  expected. 

Thus  one  may  conclude  here  that  for  the  3  KM  case,  for  the  high  sun 
cases  with  low  albedo,  the  AWS  data  agrees  well  with  the  Monte  Carlo  data 
with  a  trend  to  be  slightly  pessimistic.  For  the  high  albedo  case,  the 
AWS  data  predicts  significantly  higher  results  than  the  RRA  model.  For 
the  low  sun  case  the  AWS  data  is  significantly  low  for  the  low  albedo  case 
but  agrees  well  for  the  high  albedo  case  where  it  tends  to  be  slightly  more 
optimistic.  The  AWS  data  also  appears  to  be  much  less  sensitive  to  azimuth 
in  the  low  visibility  case  than  does  the  Monte  Carlo  data.  So  if  one  ii,  to 
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Fig.  32.  Comparison  of  Calculation  of  Contrast  Transmission  by  the  AWS 
and  RRA  Monte  Carlo  Methods.  The  Wavelength  is  .55  pm.  Solar 
Zenith  Angle  is  30°,  Albedo  is  .9.  The  Surface  Visibility  is  3  KM 
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Comparison  of  Calculation  of  Contrast  Transmission  by  the  AWS 
and  RRA  Monte  Carlo  Methods.  The  Wavelength  is  .55  ym.  Solar 
Zenith  Angle  is  30°,  Albedo  is  .9.  The  Surface  Visibility  is  3  KM 
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use  the  AWS  data  for  the  low  visibility  case,  more  care  is  needed  in  hand¬ 
ling  the  results. 

Table  VIII  shows  the  overall  performance  of  the  AWS  model  compared 
to  the  RRA  Monte  Carlo  data  for  the  two  visual  ranges.  Here  the  percentage 
of  cases  in  which  the  AWS  model  predicted  higher  values  of  contrast  trans¬ 
mission  are  shown  as  functions  of  sun  zenith  angle  and  albedo.  This  table 


TABLE  VIII 

Percentage  of  Cases  Where  AWS  Predicts  Higher  than  RRA 


Visibility 

A1 bedo 

Solar  Zenith  Angle 

0°  30°  70° 

25  KM 

.1 

83% 

77% 

79% 

.9 

90% 

96% 

89% 

3  KM 

.1 

28% 

21% 

17% 

.9 

93% 

100% 

86% 

tends  to  support  the  general  trends  noted  in  the  previous  figures.  For 
the  clear  case,  the  AWS  data  is  more  optimistic  than  the  RRA  Monte  Carlo 
data.  For  the  low  visibility  case,  AWS  is  more  pessimistic  for  the  low 
albedo  and  more  optimistic  for  the  high  albedo  case. 

Table  IX  shows  the  percentage  by  which  the  AWS  data  predicts  higher 
or  lower  than  the  Monte  Carlo  data.  From  Eq  (75)  a  negative  percentage 
indicates  that  the  AWS  data  is  predicting  higher  values  of  contrast  trans¬ 
mission  than  the  Monte  Carlo  data,  and  a  positive  percentage  indicates  AWS 
predicting  lower.  Also  shown  in  parentheses  is  the  standard  deviation  of 
the  data  set. 

From  Table  IX  it  can  be  seen  that  the  AWS  model  for  the  clear  case 
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TABLE  IX 

Percentage  Difference  Between  AWS  and  Monte  Carlo 
Predictions  of  Contrast  Transmission 


\  • 

Visibility 

1 

A1 bedo 

Solar  Zenith  Angle 

£ 

[  ? 

0° 

30° 

70° 

|  1 
\  ) 

25  KM 

.1 

-14%  (17%) 

-11%  (14%) 

-18% 

(21%) 

t  f 

1  * 

I  i 

.9 

-16%  (12%) 

-17%  (11%) 

-16% 

M 

ro 

I  * 

3  KM 

.1 

11%  (25%) 

19%  (21%) 

CO 

o 

** 

(40%) 

!  i 

.9 

-80%  (38%) 

-70%  (34%) 

00 

CVJ 

(37%) 

is  uniformly  predicting  higher  than  the  RRA  data,  and  that  the  results  are 
reasonably  comparable  for  the  six  cases.  If  one  felt  that  the  Monte  Carlo 
data  were  more  accurate,  he  might  use  the  AWS  data  because  of  its  ease  and 

I 

quickness,  and  then  scale  the  result  according  to  the  above  table.  A 
reasonable  average  might  be  -15%,  so  that  one  would  divide  the  AWS  predic¬ 
tions  by  115%  to  get  a  number  comparable  to  the  RRA  prediction. 

From  the  inconsistency  of  the  3  KM  results,  the  results  of  the  simpli¬ 
fying  assumptions  of  the  AWS  model  with  regard  to  the  number  of  scattering 
events  possible  becomes  apparent.  But  it  should  also  be  pointed  out  that 
the  results  of  the  AWS  model  are  probably  still  useful,  and  they  can  be 
corrected  from  the  above  table  to  yield  results  in  agreement  with  the  Monte 
Carlo  data. 

In  an  attempt  to  determine  the  effect  of  the  error  introduced  by 
using  different  input  altitudes  for  the  AWS  and  Monte  Carlo  data,  the  low 
altitude  point  was  eliminated  and  the  statistics  recomputed.  For  the  25  KM 
case  the  percentage  difference  became  more  negative  by  3  -  4%,  but  the 
standard  deviation  decreased  by  1  -  3%.  For  the  3  KM  case,  the  percentage 
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decreased  8  -  9%  and  standard  deviation  decreased  2-5%  for  the  low 
albedo  case.  The  percentage  decreased  11  -  17%  with  standard  deviation 
decreasing  4  -  18%  for  the  high  albedo  case. 

No  attempt  is  made  here  to  analyze  why  all  the  variation  occurs,  but 
much  is  probably  due  to  the  Monte  Carlo  data,  since  in  the  case  of  the 
solar  zenith  of  0°,  the  AWS  program  computes  the  same  contrast  transmission, 
independent  of  azimuth,  while  the  Monte  Carlo  data  does  not.  Thus  the 
standard  deviations  which  are  present  at  zenith  angle  of  0°  in  the  tables 
above  are  due  in  large  part  to  the  statistical  nature  of  the  Monte  Carlo 
data.  This  statistical  variation  might  be  expected  to  carry  over  to  the 
other  solar  zeniths  as  well . 


Comparison  to  Flight  Data 

Calculations  of  contrast  transmission  based  on  the  AWS  model  are 
compared  to  data  collected  in  flight  during  Project  Haven  View  (Ref  9). 

The  flight  measurement  techniques  were  previously  described  in  Chapter  II. 
The  processed  data  included  R* v Z , 0 , 4> ) ,  T(Z,0),  and  bRo(Zt,6,<J>) .  The  data 
was  presented  for  a  narrow  blue  and  a  narrow  red  filter,  and  the  photopic, 
and  S-20  response.  Here,  only  the  photopic  data  is  used.  With  R*,  and 
bRfl,  Eq  (58)  may  be  used  to  calculate  x  (Z,9,<j>).  The  general  weather 
conditions  and  solar  zenith  angle  are  shown  in  Table  II  in  Chapter  II. 

The  AWS  data  was  calculated  for  a  wavelength  of  .55  ym.  The  assump¬ 
tion  that  this  is  representative  of  the  photopic  region  is  made.  Several 
visibility  models  were  used  to  try  to  match  the  results  of  the  model  and 
flight  data.  The  AWS  model  with  a  surface  visibility  of  15  KM  gave  the 
best  fit  in  the  three  flights  compared  here. 

The  percentage  difference  between  the  flight  data  and  the  AWS  calcu¬ 
lation  of  contrast  transmission  is  presented  in  Table  X  for  three  flights. 
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The  standard  deviation  of  the  data  is  also  reported  in  parentheses.  Each 
entry  is  averaged  over  five  altitudes  and  seven  receiver  zenith  angles  of 
93°,  95°,  97°,  100°,  120°,  150°,  and  180°.  Slant  ranges  greater  than  15  KM 
were  ignored  in  the  calculations. 

TABLE  X 

Percentage  Difference  Between  Flight  and  AWS 
Calculation  of  Contrast  Transmission 


Flight 

(see  Chap  II) 

0° 

Azimuth 

90° 

From  Sun 

180° 

Averaged 

c-142 

-21%  (26%) 

-2.4%  (14%) 

-6%  (23.5%) 

-10%  (23%) 

c-139 

-6.5%  (16%) 

-3.8%  (17%) 

-.5%  (16%) 

-3.6%  (16%) 

c-137 

-10%  (23%) 

-6%  (24%) 

-.5%  (23%) 

-5%  (16%) 

One  may  note  first  that  the  comparison  is  reasonable  between  the  two 
sets  of  data.  The  azimuth  of  0°  is  the  largest  difference,  indicating 
that  the  AWS  model  may  not  be  handling  the  multiple  scattering  well  enough. 
Of  course  the  comparison  should  not  be  carried  too  far,  since  the  flight 
data  was  taken  under  skies  which  had  varying  amounts  of  clouds  which  the 
AWS  model  does  not  take  into  account.  Also  there  is  a  rather  wide  disparity 
between  the  reported  visibility  at  the  time  of  the  flight  and  the  visibility 
used  in  the  model  (15  KM). 

The  scattering  coefficient  data  measured  in-flight  shows  that  the 
flights  had  rather  similar  visibilities  at  the  lowest  altitude  (between 
22  and  26  KM)  but  that  a'  strong  inversion  was  present  in  all  flights  at 
about  1.5  KM,  causing  the  extinction  coefficient  to  be  nearly  constant  up 
to  that  altitude.  The  exponential  decrease  in  aerosol  extinction  coeffi¬ 
cient  up  to  4  KM  in  the  AWS  model  would  not  approximate  the  true  profile 
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closely.  The  choice  of  a  lower  surface  visibility  in  the  AWS  model  caused 
the  air-transmittance  to  match  better  so  that  the  above  reasonable  com¬ 
parison  could  be  achieved. 

The  mismatch  between  atmospheric  profiles  here  points  out  the  need  for 
continued  flight  testing  to  measure  actual  atmospheres  so  that  adequate 
atmospheric  models  can  be  constructed  for  particular  areas  and  particular 
weather  patterns.  The  AWS  program  and  most  other  atmospheric  simulation 
programs  contain  the  flexibility  to  enter  any  atmospheric  extinction  coef¬ 
ficient  profile  if  it  can  be  estimated.  If  not,  the  best  model  must  be 
used. 
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V.  Calculated  Directional  Path  Reflectance 
and  Contrast  Transmission 

In  this  chapter  of  the  report,  calculations  are  made  for  various 
combinations  of  conditions  of  sun  position  and  geometry  at  X  =  .55  pm. 

The  AWS  model  is  used  to  generate  spectral  directional  path  reflectance 
R*(Z,8,(}>)  for  each  combination  of  solar  zenith  angle,  geometry,  and  al¬ 
bedo.  This  data  is  then  displayed  in  two  forms.  Directional  path 
reflectance  is  plotted  against  azimuth  for  a  constant  altitude,  so  that 
the  effect  of  the  azimuth  of  the  path  of  sight  with  respect  to  the  sun 
may  be  determined.  Then  contours  of  constant  directional  path  reflec¬ 
tance  are  plotted  as  a  function  of  altitude  and  horizontal  range.  This 
presentation  allows  the  effects  of  receiver  zenith  angle  (or  dive  angle) 
to  be  determined. 

The  contours  of  constant  R*  are  determined  by  generating  a  table  of 
values  of  R*(Z,6,<{>)  for  10  values  of  altitude  between  1  and  10  KFT;  for 
10  values  of  0  between  95°  and  180°;  and  for  9  values  of  q>  between  0°  and 
180°.  Then  G  and  <J>  are  held  constant,  and  Z  is  determined  such  that  R* 
has  one  of  several  fixed  values.  The  interpolation  is  done  using  the 
log10R*  and  log10Z  to  improve  the  accuracy,  and  a  straight  line  interpo¬ 
lation  routine  is  used.  Then  slant  range  is  calculated  from  Range  (R*  = 
constant)  =  zjsec  0].  This  computed  range  is  plotted  against  receiver  zenith 
angle  ’nd  the  values  of  slant  range  for  constant  R*  are  connected  with  straight 
1 i nes . 

Contours  of  constant  contrast  transmission  are  generated  from  the 
R*(Z,8,<J>)  data  in  a  similar  way.  Equation  (58),  repeated  here,  is  used  to 
determine  tc(Z,6,<J>)  where  the  bR0 (Z+,0 ,4>)  data  is  taken  from  measured  data 
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in  tne  photopic  region  for  the  same  solar  zenith  angle  and  albedo. 

tc(Z,0,4>>  =  {1  +[R*(Z,6,<{))/bR0(Zt,e,({>)]}"1 

The  assumption  is  made  that  the  R*  calculated  at  .55  pm  is  represen¬ 
tative  of  a  non-spectra!  R*  for  the  entire  photopic  region.  Then  the 
contrast  transmission  is  representative  of  the  photopic  region.  The  same 
interpolation  method  as  for  R*  above  is  used  to  determine  the  contours  of 
constant  contrast  transmission. 

With  presentations  of  this  type  it  appears  that  a  good  start  may  be 
made  into  predicting  lock-on  ranges  or  detection  ranges  for  sensor  systems. 
The  presentations  should  also  be  quite  useful  for  predicting  the  areas  of 
azimuth  and  dive  angle  where  the  effects  of  non-image  forming  light  will 
be  at  a  minimum;  thus  the  seeing  conditions  should  be  the  best  for  what¬ 
ever  sensor  system  is  being  used. 

The  data  presented  here  from  the  AWS  model  calculations  is  computed 
for  the  same  solar  zenith  and  albedo  as  for  flight  data  previously  re¬ 
ported  by  Duntley  in  Ref  7,  Appendices  C  and  D.  The  DRo(Zt,0,<t>)  is  taken 
from  earlier  publications.  As  noted  in  Chapter  III  of  this  report,  other 
background  surfaces  have  also  been  measured  under  similar  conditions  and 
their  directional  reflectances  may  be  used.  Table  XI  shows  the  values  of 
the  parameters  for  which  data  is  presented.  The  figures  are  self-explanatory 
The  AWS  model  or  other  mathematical  model  easily  generates  them  for  any  situ¬ 
ation  of  interest  and  allows  parameters  such  as  visibility,  solar  zenith 
angle  and  albedo  to  be  varied  to  determine  their  effects  on  contrast 
transmission. 
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TABLE  XI 

Parameters  for  AWS  Calculations 


Solar  Zenith 
Angle 

Albedo 

Altitudes 

Visibility 

Reference  for  .  R  data 
b  o 

42° 

.06 

1-10  KFT 

23  KM 

(Ref  26:559),  also  see 
Chap.  Ill,  Figs.  16 
and  17 

78° 

.06 

1-10  KFT 

23  KM 

(Ref  23:804) 
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Contours  of  Constant  Directional  Path  Reflectance  R*  for  a  Solar  Zenith  Angle  of  42°, 
Azimuth  of  135°.  The  Albedo  is  .06  and  the  Wavelength  is  .55  yr.  Data  from  AW$  Model 


Contours  of  Constant  Contrast  Transmission  xc,  for  a  Solar  Zenith  Angle  of  42°,  and.  an  Azimuth 
of  135°.  The  Albedo  is  .06,  and  the  Wavelength  is  .55  pm.  R*  Data  is  from  AWS  Model,  t>Ro  Dat 
is  Measured  Background  Reflectance  from  Figs.  16  and  17. 


Contours  of  Constant  Directional  Path  Reflectance  R*  for  a  Solar  Zenith  Angle  of  42°  a 
Azimuth  of  180°.  The  Albedo  is  .06  and  the  Wavelength  is  .55  yjn.  Data  from  AWS  Model 


OIRECTIQNAL  PATH  REFLECTANCE 
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Fig.  49.  Contours  of  Constant  Directional  Path  Reflectance,  R*,  for  a  Solar  Zenith  Angle  of  78°,  and 
an  Azimuth  of  135°.  The  Albedo  is  .06  and  the  Wavelength  is  .55  ym.  Data  from  AWS  Model.' 
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V I ,  Inherent  Contrast  and  Spectral  Considerations 

In  this  chapter  some  techniques  for  measuring  the  contrast  between 
target  a. .5  background  will  be  discussed.  Measurements  made  at  Wright- 
Patterson  AFB  using  a  telephotometer  and  a  teleradiometer  wil1  be 
presented.  Measurements  by  photographic  means  will  also  be  shown.  The 
photographic  technique  proved  to  be  quite  convenient  and  gave  comparable 
results  to  the  telephotometer.  The  photographic  technique  was  recommended 
for  use  in  the  recent  adverse  weather  tests  of  the  Maverick  missile  in 
Europe  and  provided  much  useful  data  on  the  inherent  contrast  of  military 
vehicles.  Some  sample  data  is  presented.  Attempts  to  measure  sky-ground 
ratio  by  the  technique  suggested  by  Duntley  using  a  telephotometer,  and 
photographic  means  are  presented  and  compared  with  data  computed  from 
contrasts  measured  at  the  target  and  at  a  slant  range  of  3000  feet. 


Telephotometer  Measurements 

A  Gamma  Scientific  photometer-monochromator  with  telescopic  attach¬ 
ment  was  used.  The  photometer  has  a  photopic  response.  An  aperture  of 
6  minutes  of  arc  was  used  for  all  measurements.  The  distance  from  the 


target  and  background  was  approximately  200  feet.  The  area  subtended 
with  the  measurement  aperture  was  approximately  .35  feet,  A  2-1/2  ton 
Army  olive-drab  truck  was  measured  during  seven  separate  time  periods 
between  19  and  27  October  1971.  The  time  periods  corresponded  to  prelimi¬ 
nary  adverse  weather  flight  tests  being  conducted  by  the  Maverick  Missile 
System  Program  Office.  The  position  of  the  target  and  the  photometer  is 
shown  in  Fig.  53.  Measurements  were  conducted  both  from  the  ground  (Fig.  53a) 


and  with  the  photometer  mounted  on  a  platform  28  feet  in  the  air  (Fig.  53b) 
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i 

to  approximate  an  8°  look-down  angle,  the  approximate  dive  ancle  of  the 
aircraft. 

i 

The  reflected  radiance  in  foot-1 amberts  was  read  from  the  background 
near  the  truck,  and  from  two  areas  of  the  truck:  the  canvas  cover,  and  the 
cab.  Scans  of  different  areas  of  the  truck  indicated  that  it  was  very 
uniform  in  reflected  radiance, so  this  procedure  was  thought  to  be  repre¬ 
sentative.  The  background  was  green  grass,  about  6  inches  high  and  usually 
wee  with  dew.  Contrast  was  calculated  by  averaging  the  canvas  and  cab 
reading  as  the  target  signal  and  using  Eq  (3). 

The  results  of  the  measurements  are  shown  in  Table  XII.  The  several 
sets  of  contrast  measurements  during  a  particular  time  period  are  averaged 
and  data  measured  photographically  is  also  presented. 

Photographic  Measurements 

Photographs  of  the  target-background  scene  were  taken  during  the  same 
time  periods  as  the  telephotometer  measurements.  A  Honeywell  Pentax  camera 
with  an  f/\  *4  50-mm  fecal  length,  Takumar  lens  was  used.  Plus  X  film  was 
used  with  a  Kodak  VJratten  102  filter.  The  Wrat.ten  102  filter  is  used  to 
correct  the  response  of  the  film  to  a  phctopic  or  near  photopic  response. 

The  camera  was  located  on  the  ground  as  shown  in  Fig.  53a.  The  film 
data  was  reduced  in  the  following  way:  Prior  to  film  development,  a  stan¬ 
dard  neutral  density  step  wedge  was  exposed  on  the  film  end  as  a  check  on 
the  development  process  and  for  calibration  purposes.  The  step-wedge  has 
twenty-one  density  steps.  The  film  is  then  developed  to  a  gamma  of  approxi 
mately  1.  and  the  neutral  density  step  wedge  scanned  with  a  densitometer 
to  provide  a  D-log  E curve,  where  D  is  the  density,  and  log  E is  the  relative 
log  of  the  exposure.  A  sample  curve  is  shown  in  Fig.  54.  Although  one  can 
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actually  determine  absolute  log  E  from  the  source  used  for  the  calibration 
exposure,  this  is  not  necessary  for  the  computation  of  contrast,  if  one 
wants  the  absolute  radiance  or  luminance  of  the  target  or  background,  the 
absolute  log  E  must  be  used.  An  authoritative  treatment  of  the  photographic 
process  is  given  by  Mees  (Ref  27 : Chaps .  5,  19,  20). 

The  negative  of  the  target-background  scene  is  measured  to  yield  the 
density  of  the  target  and  background.  The  truck  image  size  was  about  5  mm. 
The  aperture  on  the  Macbeth  densitometer  was  .5  mm.  As  with  the  photometer, 
an  area  of  the  canvas  and  an  area  of  the  cab  of  the  truck  were  measured. 


Also  an  area  of  the  background  near  the  truck  was  measured.  The  D-log  E 
curve  was  then  used  to  convert  the  neutral  density  to  relative  log  expo¬ 
sure  values.  After  taking  anti -logs  to  yield  relative  exposure,  the  canvas 
and  cab  readings  from  the  truck  were  averaged  and  contrast  computed,  using 
Eq  (9). 

Averr  2d  values  of  inherent  contrast  are  presented  in  Table  XII  with 
the  telephotometer  readings  for  comparison.  In  addition,  the  sun  position 
and  sky  conditions  under  which  the  measurements  were  taken  are  listed.  The 
su  Position  data  was  computed  from  a  subprogram  of  the  Aerial  Photographic 
Energy  Model  (Ref  15).  The  visibility  is  the  estimated  visibility  in  the 
direction  of  measurement.  The  number  of  data  points  averaged  is  also  shown, 
i 

\  Radiometric  Measurements 

The  Gamna  Scientific  photometer  had  a  monochromator  attachment  which 
allowed  measurement  over  discrete  wavelength  intervals  of  width  approximately 

O 

20  A.  Measurements  were  made  from  the  platform  as  shown  in  Fig.  53b  of  the 
truck  and  the  background.  The  wavelength  intervals  covered  were  from  .45  ym 
to  .75  yn.  The  canvas  and  cab  were  measured  with  the  same  aperture  used  for  the 
telephotometer  measurements.  The  results  for  one  set  of  measurements  are 
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•*■►  shown  in  Fig.  55.  The  reflected  spectral  radiances  are  presented  in  arbi- 

■Cr 

trary  units  since  no  calibrated  source  of  known  radiance  was  available. 
However,  for  contrast  measurements  only  relative  values  are  required. 

The  data  has  been  corrected  for  the  response  of  the  S-20  phototube  which 
was  used  in  the  system  based  on  manufacturer-supplied  data.  The  contrast 
can  then  be  computed  from  Eq  (9)  where  the  relative  target  and  background 
reflected  radiances  are  computed  as  follows: 

/X=.75pm 

Nc(Z+,X,8,<J>)S(A)dX  (76) 

A=.45ym 

where  N0 (Z+,6 ,4>)  is  the  inherent  relative  reflected  radiance  in  the  range 
.45  -  .75  urn;  N0(Z+,X,6,<j>)  is  the  spectral  reflected  radiance  as  a  function 
of  wavelength  X;  and  S(A)  is  the  response  of  the  system  "seeing"  the  radia¬ 
tion  as  a  function  of  wavelength,  X. 

If  the  relative  response  of  the  eye  is  used  for  S(A),  then  the  inte¬ 
gration  carried  out,  the  result  is  the  relative  reflected  radiance  available 
to  the  eye.  The  photopic  contrast  can  then  be  computed.  This  was  done 
using  Simpson's  Rule  integration  and  data  for  three  sets  of  measurements. 

The  results  are  shown  in  Table  XII  for  comparison  with  the  photometer 
and  photographic  data.  Figure  56  shows  the  relative  radiances  available 
in  the  photopic  region.  The  response  function  for  the  eye  was  taken  from 
Brown  (Ref  28:223). 

It  can  be  noted  from  Table  XII  that  all  contrasts  were  negative  indi¬ 
cating  that  in  all  cases  the  background  was  brighter  than  the  target. 

This  is  as  anticipated  since  grass  is  normally  expected  to  have  a  reflec¬ 
tance  of  10  to  15%,  and  olive-drab  canvas  or  paint  about  5  to  10%.  The 
agreement  in  the  photometric  measurements  between  the  ground  and  platform 
measurements  is  good,  indicating  that  the  inherent  contrast  does  not  change 
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Fiq.  56.  Measured  Spectral  Radiance  for  Truck  Against  Grass 
Background.  Data  Corrected  for  Photopic  Response 
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sharply  with  small  angular  changes.  Thus  contrast  taken  from  the  ground' 
may  be  used  reasonably  for  low  dive  angles  as  long  as  pronounced  shadowing 
is  not  occuring.  Comparison  between  the  photometer  and  photographic  data 
shows  the  latter  to  be  lower  but  comparable  to  the  photometer  data, sup¬ 
porting  the  idea  that  photographs  may  be  reasonably  used  in  the  field 
environment  to  yield  inherent  contrast.  The  technique  used  here  was 
rather  simple  since  only  two  density  readings  were  made  on  the  target  and 
one  on  the  background.  This  is  justified  if  the  target  and  background 
are  relatively  uniform  as  was  the  case  here.  However,  more  complex  tar¬ 
gets  or  backgrounds  should  be  carefully  evaluated  using  a  microdensitometer 
to  determine  whether  the  simple  technique  used  above  or  a  microdensitometer 
should  be  used.  Figure  57  shows  two  scans  made  by  a  Joyce-Loebel  micro¬ 
densitometer  from  one  of  the  frames. •  An  aperture  „ .  25  pm  by  615  ym  was 
used.  The  large  density  change  between  target  and  grass  background  can 
be  clearly  seen.  But  note  also  the  relatively  smooth  density  readings 
from  the  canvas  near  the  top  of  the  truck  (Figure  57a),  but  the  more 
irregular  readings  near  the  bottom  (Figure  57b). 

Figures  58  and  59  show  microdensitometer  maps  of  a  2-1/2  ton  truck 
and  a  Chapparal  weapons  carrier.  They  were  generated  by  computer  from 
negatives  using  the  same  camera  and  film  as  above.  These  were  taken  in 
Germany  during  the  recent  adverse  weather  tests  of  the  Maverick  missile. 

Only  the  target  vehicle  and  the  surface  background  are  displayed,  the 
relative  density  values  being  indicated  by  the  darkness  or  lightness  of 
the  areas.  Note  the  uniformity  of  the  truck  with  certain  darker  sha¬ 
dows  near  the  wheels,  while  the  chapparal  weapons  carrier  is  much  less 
distinct  and  blends  with  the  background.  Note  also  that  this  is  a  "positive" 
while  the  densities  shown  in  Figs.  57a  and  57b  are  "negatives."  To  get 
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Fig.  59.  Computer  Generated  Picture  of  Chapparal  Weapons  Carrier  Used  D’  •  >v-  ck  Adverse  Weather 

Testing  in  Germany.  (Picture  was  generated  by  the  same  techn  58.  Note  the  lower 

contrast  presented  by  this  target  in  comparison  to  the  2-1/2  .  ui'ts  also  the  very 

bright  strip  and  the  dark  areas  between  the  treads.) 
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meaningful  inherent  contrasts  on  the  Chaparral  it  might  be  necessary  to 
go  to  the  microdensitometer  technique.  Inherent  contrasts  achieved  in 
Europe  during  the  Maverick  tests  will  be  reported  in  a  Flight  Test  report 
to  be  published  by  the  Maverick  System  Program  Office. 

Spectral  Considerations 

Throughout  the  report,  the  fact  that  the  radiometric  and  atmospheric 
quantities  are  spectral  in  nature  has  been  mentioned,  brf  "hing  has  been 
done  further  than  to  assume  that  a  wavelength  of  .55  ytr  'esentative 

of  the  photopic  region  and  state  that  this  is  the  area  o  .ern.  For 
an  electro-optical  sensor,  this  is  not  always  the  case..  If,  for  example, 
the  measured  spectral  radiances  from  the  truck  and  background  presented  in 
Fig.  55  are  integrated  over  wavelength  using  the  spectral  response  of  the 
vidicon.  Fig.  61  results.  The  response  function  for  the  vidicon  was  taken 
from  a  Hughes  technical  report  with  the  vidicon  corrected  by  a  yellow  haze 
filter.  The  response  function  of  the  vidicon  is  shown  in  Fig.  60.  From 
Fig.  61  it  can  be  seen  that  the  vidicon  sees  more  contrast  in  the  red  end 
of  the  spectrum,  having  sensitivity  into  the  near  IR-  he  inherent  con¬ 
trast  was  computed  for  the  three  data  sets  irom  Table  MI  where  radiometric 
measurements  were  available.  The  integration  was  carried  out  from  .45  ym 
to  .75  ym  using  the  measured  data,  and  the  data  at  the  end  points  was  then 
extrapolated  to  .4  ym  and  .8  ym  as  shown  in  Fig.  55.  The  integration  was 
then  carried  out  from  .4  ym  to  .8  ym.  The  results  are  presented  in  Table 
XIII.  Note  from  the  table  that  more  contrast  exists  for  the  vidicon  than 
for  the  eye  and  that  integration  to  longer  wavelengths  is  helpful  for  im¬ 
proving  the  contrast  available.  From  the  response  curve  for  the  vidicon, 
integration  to  about  .85  ym  would  be  justified  and  would  probably  show  that 
even  more  inherent  contrast  is  available  to  the  vidicon  than  to  the  eye 
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Fig.  60.  Relative  Spectral  Response  of  Maverick  Vidicon 
And  With  Yellow  Haze  Filter 


A.  ^ 


RADIANCE  (ARBITRARY  UNITS) 


Fig.  61.  Measured  Spectral  Radiance  for  Truck  Against  Grass 
Background.  Data  Corrected  for  Vidicon  Response. 
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TABLE  XIII 


Inherent  Contrast  Available  t.o  the  Eye  and  to  a  Vidicon 


Date, 

Time  (EOT) 

Inherent  Contrast 

Photopic  Response  Vidicon  With 

.45-. 75  ym  .4-. 8  ym  .45-. 75  ym 

Filter 
.4-. 8  ym 

26  Oct  71 
1130-1200 

-.66 

-.56 

-.68 

-.70 

26  Oct  71 
1455-1530 

-.66 

-.66 

-.69 

-.71 

27  Oct  71 
1020-1050 

-.89 

■  .89 

-.90 

-.91 

for  this  target  background  situation.  In  ;his  situation  advantage  exists 
for  this  vidicon  over  the  eye  with  regard  to  inherent,  contrast  available. 

The  atmospheric  contrast  transmission  is  also  more  favorable  for  the 
longer  visible  wavelengths.  The  RRA  Monte  Carlo  data  for  an  isotropic 
source  with  albedo  of  .1,  receiver  at  .5  KM,  and  a  receiver  zenith  angle 
of  98°  is  given  in  Table  XIV.  This  data  would  indicate  a  higher  apparent 
contrast  available  at  the  sensor  for  the  vidicon  than  for  the  eye,  due  to 
increased  contrast  transmission.  Thus  in  any  carefu.  analysis  in  which 
predictions  of  sensor  performance  are  to  be  made,  the  spectral  considera¬ 
tions  should  be  taken  into  account.  Then  the  apparent  contrast  should  be 
determined  from 

cr(z,0,4>)  =  C  Tcu,z,e,4>)c0u,z+/e,<f>)dX  (77) 

A  of 

sensor 
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TABLE  XIV 

Spectral  Contrast  Transmission  From  Monto  Carlo  Data 


Wavelength 

.45  ym 
.55  ym 
.65  ym 
.75  ym 
.85  ym 


Contrast  Transmission 


and  the  required  apparent  contrast  of  the  system  should  be  specified  in 
terms  of  the  response  if  the  system. 

Apparent  Contrast  Measurements 

During  the  periods  while  measurements  were  made  at  Wright-Patterson 
AFB,  a  camera  with  a  1000  mm  focal  length,  //12  lens  v/as  positioned  approxi¬ 
mately  3000  feet  from  the  truck.  The  camera  was  located  along  the  same 
azimuth  as  the  phonometer  and  inherent  contrast  camera „  It  was  on  a  small 
hill  giving  a  look-down  angle  of  .8°.  The  camera  also  used  Kodak  Plus-X 
film  and  a  Wratten  102  filter.  Photos  of  the  truck-background  scene  were 
taken  simultaneously  with  those  taken  at  close-range  and  the  film  processed 
and  analyzed  in  the  same  way  to  yield  contrast.  Thus  the  apparent  contrast 
at  this  position  was  determined.  The  contrast  transmission  could  then  be 
calculated  from  the  ratio.  The  results  are  shown  in  Table  XV.  From  the 
estimated  visibility  and  the  measured  contrast  transmission  the  sky-ground 
ratio  was  calculated  using  Eq  (43).  These  values  are  shown  also  in  Table  XV. 
In  an  attempt  to  determine  whether  the  sky-ground  ratio  could  be  measured 
in  the  way  indicated  by  Duntley,  photometer  measurements  of  the  sky  and 
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Numbers  in  parentheses  are  the  number  of  samples  averaged 
H  indicates  haze  obstructing  visibility 

Negative  sign  on  contrast  indicates  target  darker  than  background 
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the  ground  were  taken  in  the  appropriate  direction.  In  this  case  as  can 
be  seen  from  Eqs  (37)  and  (38),  the  situation  reduces  to  the  case  where 
the  sky  should  be  measured  in  the  same  direction  as  the  sensors  are 
looking  at  the  truck.  The  averaged  sky-ground  ratios  by  this  technique 
are  shown  in  Table  XV  with  the  number  of  samples  averaged.  Here  as  in 
all  cases  the  "ground"  value  is  the  reflected  radiance  from  the  direc¬ 
tion  of  the  target  along  the  inclined  path  to  the  sensor.  Similar 
measurements  were  made  from  the  inherent  contrast  photos  of  sky  and 
ground  in  the  direction  of  the  target.  The  averaged  values  are  shown  in 
Table  XV. 

From  this  tar.le  it  can  be  seen  that  the  measured  contrast  trans¬ 
missions  and  visibilities  yield  significantly  lower  values  of  sky-ground 
ratio  except  in  the  case  where  the  sky  is  clear.  Thus  the  results  indi¬ 
cate  that  when  there  is  an  overcast  or  broken  sky  condition,  the  sky-ground 
ratio  cannot  be  measured  directly.  With  the  clear  sky  condition,  signifi¬ 
cant  haze  was  present,  causing  the  sky-ground  ratio  values  to  differ  from 
what  could  be  measured  by  the  technique  suggested  by  Duntley.  Two  values 
are  shown  for  the  values  of  sky-ground  ratio  computed  from  inherent  and 
apparent  contrast  due  to  the  uncertainty  in  the  estimated  visibility.  The 
larger  value  is  based  on  studies  which  show  that  the  meteorologist  typi¬ 
cally  estimates  about  75%  of  the  true  meteorological  range.  True  meteoro¬ 
logical  range  is  defined  as  that  range  where  a  large  black  object  is  just 
visible  against  the  horizon  sky.  The  apparent  contrast  threshold  of  the 
eye  is  assumed  to  be  .02. 
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VII.  Conclusions  and  Recommendations 


As  a  result  of  the  study  several  conclusions  may  be  drawn: 

1.  The  prediction  of  contrast  transmission  by  using  the  sky-ground 
ratio  and  optical  standard  atmosphere  must  be  modified  to  take 
into  account  the  more  realistic  values  of  sky-ground  ratio  and 
atmospheric  air-transmittance  given  in  Chapter  II. 

2.  The  directional  path  reflectance  R*,  proposed  more  recently  by 
Duntley,  is  a  much  better  single  parameter  for  use  in  evaluating 
the  condition  of  the  atmosphere  for  seeing.  It  is  representative 
of  the  solar  zenith, viewing  angle,  atmosphere  and  general  scene 
albedo.  R*  may  be  used  to  evaluate  the  best  azimuths  and  dive 
angles  for  any  particular  task.  The  value  may  be  obtained 

from  flight  measurements  or  from  model  calculations  in  which 
assumptions  are  made. 

3.  Comparison  of  the  RRA  Monte  Carlo  model  and  the  AWS  model  for 
the  prediction  of  contrast  transmission  showed  that  the  AWS  model 
predicted  higher  contrast  transmission  for  the  25  KM  case  by  15%. 
For  the  3  KM  case,  the  AWS  data  predicted  lower  contrast  trans¬ 
mission  for  the  low  albedo  case  but  mugh  higher  contrast  trans¬ 
mission  tor  the  high  albedo  case.  This  inconsistency  is  probably 
explained  by  the  fact  that  the  AWS  model  accounts  for  only  first 
and  second  order  scattering,  while  the  atmospheric  optical  path 
becomes  very  long  for  low  sun  angles  and  low  visibilities.  The 
Monte  Carlo  data  takes  into  account  higher  order  scattering. 

The  AWS  model  or  other  models  of  its  type  are  optimized  for  fast 
running  on  the  computer.  They  provide  a  valuable  tool  for 
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generating  the  needed  data  for  making  predictions  of  atmospheric 
effects  on  contrast  transmission. 

4.  The  inherent  directional  background  reflectance  was  shown  to  be 
very  important  in  the  prediction  of  contrast  transmission.  A 
general  tendency  observed  in  measured  data  for  the  directional 
reflectance  to  increase  at  the  lower  dive  angles  tends  to  par¬ 
tially  compensate  for  a  decrease  in  air- transmittance.  There  Is 
also  generally  an  increase  near  the  specular  and  anti -solar  angles. 
Thus  the  contours  of  constant  contrast  transmission  may  be  signi¬ 
ficantly  modified  when  compared  to  those  resulting  from  the 
assumption  of  a  purely  Lambertian  surface. 

5.  Techniques  were  investigated  for  the  measurement  of  contrast. 

The  photographic  technique  was  found  to  be  comparable  to  the 
photometer  technique.  Simple  densitometer  techniques  proved 
adequate  for  relatively  uniform  target  and  background  scenes. 

This  photo  technique  provides  a  very  simple  method  of  making  in¬ 
herent  or  apparent  contrast  measurements  during  flight  tests. 

6.  Attempts  to  measure  the  sky-ground  ratio  by  the  technique 
suggested  by  Duntley  in  1948  proved  to  be  a  failure  when  the 
skies  were  broken  to  overcast.  The  agreement  was  closer  when 
the  skies  were  clear  although  moderate  haze  was  present.  This 
confirmed  Duntley 's  observation  that  the  equilibrium  radiance  was 
equal  to  the  radiance  along  the  horizon  only  in  the  case  of 
clear  skies. 

It  is  recommended  that  several  actions  be  taken  to  improve  predictive 
techniques  and  to  validate  the  models  through  testing: 

1.  The  AWS  model  and  the  RRA  Monte  Carlo  model  should  be  compared  to 


135 


,-t,^._..^1 _ ^  n . .  ,,. 

GEP/PH/72-4 


<  ► 

v 


other  of  the  existing  mathematical  models  and  to  flight  data 
to  validate  the  models  and  determine  their  limitations.  A  data 
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set  has  been  requested  from  the  University  of  Michigan.  This 
will  allow  comparison  under  the  same  conditions  as  described 
in  Chapter  IV. 

A  desirable  outcome  of  this  research  would  be  to  generate  a 
closed  form  analytic  model  which  predicts  contrast  transmission 
based  on  the  geophysical  and  weather  parameters.  The  computer 
models,  particularly  the  fast  models,  provide  a  good  opportunity 
to  generate  large  amounts  of  data  under  varying  parameters  for 
use  in  multiple  regression  analysis.  Thus  the  effects  and  mag¬ 
nitudes  of  the  effects  of  different  parameters  can  be  investigated 
and  perhaps  an  analytic  model  evolved. 

Models  are  only  as  good  as  the  correct  predictions  they  make. 

They  must  be  checked  against  measured  data.  There  is  a  clear 
lack  of  reliable  meteorological  and  optical  data  during  almost 
all  flight  tests  of  electro-optical  systems.  This  is  largely 
because  the  requirements  of  the  atmospheric  physicists  and  the 
system  program  engineers  do  not  overlap.  The  atmospheric  data 
is  usually  an  afterthought,  with  test  engineers  accepting  the 
standard  meteorological  information  available.  It  is  only 
through  close  planning  at  high  levels  that  the  two  requirements 
can  be  overlapped.  This  requires  that  flight  tests  be  planned 
far  enough  in  advance  that  the  atmospheric  aspects  may  be  con¬ 
sidered  and  integrated  into  the  program. 

Many  of  the  effects  important  to  atmospheric  contrast  transmission 
can  be  studied  without  the  particular  electro-optical  sensor  being 
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present.  For  example,  the  decoupling  of  the  target  from  the 
atmospheric  contrast  transmission  implies  that  inherent  contrasts 
may  be  measured  for  typical  military  targets  against  typical  back¬ 
grounds  without  the  particular  aircraft  and  sensor  flying  around. 
This  means  that  much  more  data  can  be  achieved  on  inherent  con¬ 
trast  under  a  wide  variety  of  conditions,  and  at  an  enormously 
decreased  price.  Measurements  like  these  would  also  provide  In¬ 
herent  directional  background  reflectances  which  are  important 
to  contrast  transmission.  This  requires  that  continued  emphasis 
be  given  to  applied  research  programs,  so  that  the  answers  will 
be  available  for  future  systems  in  the  concept  stages  and  not 
later  in  the  testing  phase. 
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Appendix  A 

Definitions  of  Contrast 

Three  definitions  for  spectral  contrast  are  examined  here  along  with 
related  expressions  for  spectral  contrast  transmission.  All  are  based  on 
the  equation  relating  the  spectral  radiance  at  an  altitude  Z,  along  a 
slant  path  with  receiver  zenith  angle  9  and  azimuth  from  the  sun  <}>  to  the 
spectral  radiance  leaving  the  source  in  that  direction. 

N  (Z,0,<|O  =  N  <Z.,e,4>)T(Z,9)  +  N*(Z,9,<J>)  (A.l) 

r  o  “ 

Universal  Contrast 

The  definition  of  universal  contrast  was  given  in  Chapter  II,  Eqs 
(8)  and  (9),  and  for  universal  contrast  transmission  in  Eqs (13)  and  (58). 
Using  Eq  (10), 


Cr(Z,0,<j>)  =  Co(Z,0,d')Tc(Z,0,<i>) 
using  Eqs  (8)  and  (13) 


(A. 2) 


cr(z,e,«j>)  = 


+No(Z+,0,cf>)  -  bNo(Z+,0,<j» 


bW6^ 


bN0(Z+,9,<j)}  +  N*(Z,0,<j>)/T(Z,9) 


(A-3) 


+No(Z+>0,({))  -  bN0(Z+,6,<{>) 


cr(z,9,4))  =  N#(z  e^)  +  N*(z,e,<}>)/T(zf0) 


(A-4) 


dividing  each  term  in  the  numerator  and  denominator  by  H(Zt,d)/ir ,  the 
spectral  irradiance  on  an  upward  facing  horizontal  surface  at  the  target 
altitude,  and  recognizing  that  from  Eq  (58)  these  quantities  are  bR<>  (Z+,0,<J>), 
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+R„ (Zt,6,4>)  and  R*<Z,e,4>) ,  the  directional  spectral  reflectances  of  back¬ 
ground,  target,  and  path  respectively,  then 


+R0<z,e,<f>)  -  bR0(z,e,<j>) 
cr(z,e,4>)  -  bRo(Z>9><{))  +  R*(z,e,<})) 


(A. 5) 


Thus  a  useful  relation  between  apparent  spectral  contrast  and  R*, 
bR0  and  +Ro  is  obtained. 

Contrast  Modulation 

Inherent  spectral  contrast  modulation  is  defined  as 


+N0 (Z,e,4>)  -  bN0CZ,8,4>) 

+n0  (z,e,4>)  +  bN0(z,e,<j>) 


(A. 6) 


Developing  an  expression  for  spectral  contrast  modulation  at  altitude 
Z,  receiver  zenith  angle  0  and  azimuth  <j>,  we  get 


+N0(Z,e,<j>)T(Z,6)  -  bN0(Z,e,<j>)T(Z,( 


-j-  [)  • WWj  —  /  ^  g r  r  »  - r  -  - 

cr^(z,o,(j))  =  +No(Ztj9  4. ) tTz , e , <t> )  +  bN0(z+,e,4>)  +  2N*(z,e,d>)  {A'7) 


dividing  by  T(Z,6) 


+N0(Zt,9,4>)  -  bNo(Z+,0,4>) 
+No(Z+,0,<j>)'+  bNQ  (Z+,6,4>)  +  2N*(Z,0,4>)/T  (Z,0) 


(A. 8) 


Tc^(Z,0,<j.)  =  +No(Z^>6>4))  +  bN0(Z+,e,4>)  +  2N*(Z,0,<j>)/T(Z,O>  {A*9) 


+No(Z+,0,4>)  +  bNo(Z+,0,4>) 


dividing  each  term  in  the  denominator  by  H(Z,d)/ir  we  get 
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T  (Z,0,4>) 
cm 


_ 1 _ 

,  .  _ R*(z,et(j>) 

[fR  (zt,e,4>)  +  bRo(Z+,0,4>):/2 


(A. 10) 


then 


xcjz,e,(j))  -  1  +  R*(Z>e,<j,)/R(Z+,e,<j>) 


(A.n) 


where  Ro(Zt,0,<|>)  is  the  average  spectral  directional  reflectance  of  the 
target  and  background. 

Then 


Cr  (Z,e,<j>)  =  Tc  (Z,0,<{))CO  (z+,e,4>)  (a. 12) 

m  m  m 


_  AWiVf  (A.13) 

rm  1  +  R^<Z,e,4>)/R0  (z+,0,d>) 

Contrast  Ratio 

Spectral  contrast  ratio  is  defined  as  the  ratio  of  target  spectral 
radiance  to  background  spectral  radiance. 

Then 

,No(Z+,0,4>) 

C0  (Z+,0,<J>)  =  n  (Z  0,$)  ^A*14) 

cn  b  0  T 


and 


Cr  (Z,0,(j>) 
CR 


+No(Z+,0,(j))T(Z,e)  +  N*(Z,0,<t>) 
bNo(Z+,0,4>)T(Z,0)  +  N*(Z,0,d>) 


dividing  by  T(Z,0)  and  chen  by  H(Z,d)/ir, 


(A. 15) 
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+R0(zt,e,«j>)  +  R*(zfe.<j>) 
cr^(z+,e,d))  =  bRo(Z+>e,4))  +  r*<z,m> 


and 


(A. 16) 


tc  (z,e,4>) 


i  +  f^(z,e,4>)/+R0(z,e,<{>) 
1  +  R*(Z,9,4>)/bRo(Z,e,4>) 


(A. 17) 


We  note  that  useful  relations  have  been  developed  between  three 
Jefinitions  of  contrast  and  the  directional  path  reflectance  R*(Z,6,<J>), 
and  the  inherent  spectral  directional  reflectances  of  target  and  back¬ 
ground.  From  the  flight  measurement  point  of  view,  the  universal  contrast 
definition  is  the  most  useful, since  we  see  that  only  in  this  definition  is 
the  spectral  contrast  transmission  independent  of  the  target. 
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Appendix  B 

Computation  of  Scattering  Angle 
To  compute  angle  £  between  two  line',  in  space, 

cos  0.  =  cos  otj  cos  a2  +  cos  0j  cos  02  +  cos  yx  cos  y2  (B.l 


where 


cos 

al 

_  Xl 

d 

sin  0O 

cos  <J)0 

d 

d 

cos 

a 

l 

=  sin 

9o 

cos  4>0 

y  2 

d  ; 

sin  0O 

sin  <*><, 

cos 

=  d”  = 

d 

cos 

6i 

=  sin 

6o 

sin  <j)0 

z 

d 

cos  0- 

cos 

1-° 

11 

0 

d 

=  cos  0O 

then 


cos  ax  =  sin  0O  cos  <f>0 
cos  Bj  =  sin  0O  sin  <J>0 
cos  Yi  =  cos  0O 


(B.2 

(B.3 

(B.4 


For  an  observer  at  receiver  zenith  angle  0  and  azimuth  <£_, 


T2 

cos  a 

2 


180  -  0 

x2  d2  sin( 180  -  0)  cos  £ 


cos  a  =  s  i  n  ( 1 80  -  0)  cos  4> 


t  irft «“**~ - -  *■*• /&■*.! 
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in  terms  of  angles  from  the  sun  in  the  yz  plane 

<j>  =  <j)  -  90°  =  -  (90  -  4>) 


than 


cos  B  =  -  sin  90  sin  0  cos  <j>  -  cos  0 


o  cos 


The  scattering  angle  0  is  180°  -  8,  so 
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Appendix  C 

Calculated  Sky-Ground  Ratios 
From  RRA  Monte  Carlo  Data 

This  appendix  presents  sky-ground  ratios  calculated  from  the  RRA 
Monte  Carlo  data.  The  sky-ground  ratio  is  calculated  from  Eq  (43)  by 
solving  for  K  in  terms  of  tc  and  the  air-transmittance  T  =  exp(  -aoR). 

Both  the  contrast  transmission  and  the  air-transmittance  are  known  from 
Ref  8.  Due  to  the  variations  of  K  with  altitude,  all  values  for  a  parti¬ 
cular  receiver  zenith  angle  have  been  averaged  over  the  altitudes  for 
which  data  was  computed  (.5  KM  to  10  KM).  For  a  solar  zenith  angle  of 
G°,  only  the  azimuth  averaged  values  of  K  are  presented,  since  for  this 
solar  zenith,  any  variation  with  azimuth  is  meaningless  and  due  to  the 
statistical  variation  of  the  Monte  Carlo  method. 

Values  of  sky-ground  ratio  for  slant  paths  in  excess  of  12  KM  for 
the  3  KM  surface  visibility  case  and  50  KM  for  the  25  KM  surface  visibi¬ 
lity  case  have  been  excluded  from  the  altitude  averaging  This  is  done 
because  an  increase  in  error  in  the  calculation  of  air-transmittance  for 
long  slant  paths  causes  a  sharp  increase  in  the  computed  K  value.  The 
exclusion  is  justified  by  the  rather  academic  usefulness  of  values  at 
slant  ranges  greater  than  12  or  50  KM  for  the  particular  visibilities 
quoted. 
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TABLE  C-I 

Computed  Sky-Ground  Ratios  from  RRA  Monte  Carlo 
Data  for  3  KM  Visibility 

All  entries  have  been  averaged  over  altitude.  Values  for  slant  ranges 
greater  than  12  KM  have  been  ignored. 


Receiver 

Zenith 

Angle 


Albedo 


e 

.1 

.3 

.6 

.9 

Wavelcr,..! 

171 

3.4 

1.9 

1.5 

1.4 

X  =  .55  pm 

164 

3.6 

2.0 

1.6 

1.4 

Overcast 

157 

3.5 

2.0 

1.6 

1.4 

148 

3.9 

2.1 

1.6 

1.5 

138 

4.5 

2.3 

1.7 

1.5 

130 

5.1 

2.5 

1.8 

1.6 

123 

6.6 

2.9 

2.0 

1.7 

116 

6.3 

3.0 

2.1 

1.7 

110 

11. 

4.0 

2.5 

2.0 

104 

11. 

4.4 

2.7 

2.1 

98 

3.4 

1.5 

.96 

.79 

93 

3.2 

1.3 

.77 

.61 

Wavelength 

171 

3.9 

2.3 

1.9 

1.8 

X  =  .65  ym 

164 

4.75 

2.6 

2.1 

1.9 

Overcast 
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<P  TABLE  C-Ii 

Computed  Sky-Ground  Ratios  from  RRA  Monte  Carlo 
Data  for  25  KM  Visibility 


All  entries  have  been  averaged  over  altitude.  Values  for  slant  ranges 
greater  than  50  KM  have  been  ignored. 


Receiver 

Zenith 


Albedo 


Angle 

6 

.1 

.3 

.6 

.9 

Wavelength 

171 

2.6 

1.4 

1.1 

1.0 

X  =  .55  ym 

164 

2.5 

1.3 

1.0 

.90 

Overcast 

157 

2.7 

1.4 

1.0 

.92 

148 

1.8 

1.3 

1.0 

.90 

138 

2.9 

1.4 

1.0 

.85 

130 

3.5 

1.6 

1.2 

1.0 

123 

3.8 

1.7 

1.2 

10 

116 

3.9 

1.7 

1.1 

.96 

110 

4.3 

1.8 

1.2 

.94 

104 

4.9 

1.9 

1.2 

.96 

98 

5.1 

1.9 

1.1 

.77 

93 

5.1 

2.0 

1.3 

1.0 

Wavelength 
X  =  .65  ym 

Overcast 


171 

2.3 

1.2 

.93 

.83 

164 

2.4 

1.3 

1.0 

.93 

157 

2.0 

1.0 

.76 

.70 
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TABLE  C-III 

Computed  Sky-Ground  Ratios  from  RRA  Monte  Carlo 
Data  for  3  KM  Visibility 


All  entries  have  been  averaged  over  altitude, 
greater  than  12  KM  have  been  ignored. 

Wavelength  X  =  .55  ym;  Solar  Zenith  Angle  60 

Values  for  slant  ranges 

=  0° 

Receiver  Zenith  Angle 
e 

Azimuth  Averaged 

Albedo  =  .1 

171 

2.6 

164 

2.7 

157 

3.0 

148 

2.7. 

138 

2.8 

130 

2.5 

123 

3.3 

116 

2.9 

110 

3.9 

104 

3.9 

98 

1.9 

93 

1.6 

Albedo  =  .9 

171 

1.3 

164 

1.4 

157 

1.4 

148 

1.3 

138 

1.4 

130 

**  1.4 

123 

1.5 

116 

1.3 

110 

1.6 

104 

1.4 

98 

.6 

93 

.4 
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TABLE  C-IV 

Computed  Sky-Ground  Ratios  fr'rni  RRA  Monte  Carlo 
Data  for  25  KM  Visibility 

All  entries  have  been  averaged  over  altitude.  Values  for  slant  ranges 
greater  than  50  KM  have  been  ignored. 

Wavelength  X  =  .55  ym;  Solar  Zenith  Angle  60  =  0° 


Receiver  Zenith  Angle 
6 


Azimuth  Averaged 


tapmsb***  -t-^ 


Y'l 


1 
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TABLE  C-V 


Computed  Sky-Ground  Ratios  from  RRA  Monte  Carlo 
Data  for  3  KM  Visibility 


All  entries  have  been  averaged  over  altitude.  Values  for  slant  ranges 
greater  than  12  KM  have  been  ignored. 


Wavelength  A  =  .55  ym;  Solar  Zenith  Angle  8o  =  3P° 
R?cej??r  I  Azimuth 


Receiver 

Zenith 

Angle 

6 


m 


»T^V  I?w7,)7i;^f?\^^!l».  ikis»3!\  v"v^*/,»  .*  ‘~?-».V>?>"-W»' 
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TABLE  C-VI 


Computed  Sky-Ground  Ratios  from  RRA  Monte  Carlo 
Data  for  25  KM  Visibility 


All  entries  have  been  averaged  over  altitude.  Values  for  slant  ranges 
greater  than  50  KM  have  been  ignored. 


Wavelength  X  =  .55  ym;  Sola"  Zenith  Angle  60  =  30° 
R!cel!ff  I  Azimuth 


Receiver 

Zenith 

Angle 

6 


22.5° 

67.5° 

112.5° 

157.5° 

2.6 

2.1 

3.5 

1.6 

1.8 

2.4 

1.8 

3-1 

1.8 

1.2 

1.4 

1.8 

1.7 

1.4 

1.9 

2.0 

1.5 

1.9 

1.7 

1.7 

2.0 

1.7 

1.7 

1.8 

2.2 

1.6 

2.1 

2.4 

1.6 

1.7 

1.6 

1.9 

1.9 

1.8 

2.2 

1.9 

1.4 

1.6 

1.5 

1.9 

1.7 

2.1 

2.3 

1.3 

3.4 

1.9 

1.6 

1.4 

•gw.*""!! 


Sin 


TABLE  C-VII 

Computed  Sky-Ground  Ratios  from  RRA  Monte  Carlo 
Data  for  3  KM  Visibility 

All  entries  have  been  averaged  over  altitude.  Values  for  slant  ranges 
greater  than  12  KM  have  been  ignored. 
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TABLE  C-VIII 

Computed  Sky-Ground  Ratios  from  RRA  Mont°  Carlo 
Data  for  25  KM  Visibility 


All  entries  have  been  averaged  over  altitude.  Values  for  slant  ranges 
greater  than  50  KM  have  been  ignored. 


Wavelength  X  =  .55  vim;  Solar  Zenith  'Migle  60  =  70° 


Receiver 

Zenith 

Angle 

0 

Azimuth 

22.5° 

67.5° 

112.5° 

157.5° 

Average 

Albedo  = 

.1 

171 

2.3 

4.8 

3.7 

1.8 

3.1 

164 

2.5 

2.2 

4.4 

2.8 

2.8 

157 

2.9 

4.0 

3.3 

4.3 

3.6 

148 

3.2 

2.9 

3.9 

3.6 

3.4 

138 

3.3 

3.8 

2.8 

3.8 

3.4 

130 

5.8 

3.6 

3.9 

4.1 

4.4 

123 

6.0 

3.9 

4.1 

3.8 

4.5 

116 

8.8 

5.3 

3.8 

5.4 

5.8 

110 

9.8 

4.2 

4.4 

5.6 

104 

11.5 

5.5 

3.4 

5.4 

6.5 

98 

17.0 

4.6 

3.3 

4.8 

7.4 

93 

15.0 

5.3 

7.4 

3.6 

7.9 

Albedo  = 

.9 

171 

.37 

1.51 

1.52 

.74 

1.0 

164 

.16 

1.2 

1.1 

.98 

157 

.50 

1.4 

1.2 

1.2 

148 

.90 

.73 

1.4 

.93 

) 

138 

.89 

1.1 

.77 

.86 

.90 

130 

1.5 

.89 

1.1 

.99 

1.1 

123 

1.2 

.78 

1.1 

.83 

.98 

116 

1.25 

1.04 

.88 

1.1 

1.1 

110 

1.4 

.92 

1.2 

1.0 

1.1 

104 

1.7 

.92 

.82 

1.0 

1.1 

98 

2.0 

.86 

.74 

.81 

1.1 

93 

2.0 

.89 

1.2 

.81 

1.2 
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